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Inertial  waves  excited  in  the  aixed  layer  by  hurricane 
Prederic,  had  horizontal  scales  of  approxiaately  1  to  2 
tiees  the  taroclinic  Rossfcy  radius  of  deforaation  (50  km)  of 
the  first  eode  near  the  EeSoto  Canyon.  Initially,  energy 
propagated  vertically  at  about  1.25  ka/d  and  horizontally  at 
80  ka/d.  These  waves  spun  down  over  enfolding  scales  of 
four  inertial  periods  as  energy  propagated  vertically  at  270 
a/d  and  horizontally  at  30  ka/d.  Inertio-gravity  waves  in 
the  deep  thertocline  had  horizontal  scales  of  25  to  50  ka 
and  vertical  scales  approxiaately  equal  to  the  water  depth. 
The  energy  of  these  waves  was  dcainated  by  the  barotropic 
node  with  soac  contributions  frca  aodes  1  and  2.  These 
waves  were  not  adaitted  to  the  shelf  region  because  the  bct- 
toa  slope  was  greater  than  the  slope  of  the  internal  wave 
characteristics. 

The  aean  flew  followed  the  isobaths  at  all  levels,  but 
it  was  in  the  opposite  direction  in  the  bottoa  layer.  The 
aean  flow  initially  decreased  alcng  the  eastern  boundary  of 
the  canyon  as  the  store  fencing  readjusted  the  flow.  Near* 
bottca  teaperature  variations  of  4  0c  were  associated  with 
the  store  surge  and  advection  in  the  along-track  direction, 
particularly  along  the  north  ria  of  the  canyon. 
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X.  INTROD OCTION 


"  the  doer  alone  learn eth  n 
Nietzsche . 

The  oceanic  response  tc  intense,  transient  ataospheric 
events,  such  as  hurricanes,  is  dcninated  by  robust  inertial 
wave  excitation  in  the  nixed  layer  and  the  subsequent  propa¬ 
gation  and  dispersion  of  inertio-gravity  waves  in  the  ther- 
aocline.  The  study  of  these  forced  waves  has  been  largely 
through  nuaerical  and  analytical  investigations,  because  of 
the  lack  of  sufficient  data  to  resolve  the  scales  and  ener¬ 
getics  of  the  notion.  Hence,  observations  of  ocean  currents 
and  teaperatures  daring  a  hurricane  are  neccessary  for  vali¬ 
dating  aodels  and  theories  dealing  with  forced  inertio-grav¬ 
ity  waves. 

A  conprebensive  data  set  of  ocean  observations  was  col¬ 
lected  by  the  0.  S.  Naval  Oceanographic  Office  (NAVOCEANO) 
during  the  passage  of  hurricane  Frederic  in  the  suaaer  of 
1979  (Shay  and  Tanul,  1980).  NAVOCEANO  deployed  three 
noored  taut-wire  current  inter  arrays  in  the  northern  Gulf 
of  Hexico  southeast  of  Hobile,  Alabaaa.  During  the  period 
of  deployaent,  hurricane  Frederic  passed  within  80  to  130  ka 


- 


of  the  aria;  sites  (Fig.  1)  at  2  100  GHT  12  Sept9aber.  Five 
hours  later#  Frederic  Bade  landfall  at  Dauphine  Island,  Ala- 
baaa.  During  and  subsequent  to  the  passage  of  Frederic,  the 
Haticnal  Hurricane  Research  Division  dropped  several  expend¬ 
able  bathytheracgraphs  (AXE?)  froa  reconnassiance  aircraft 
in  the  area  cf  the  current  aeter  deployaents  (Black,  1983). 


Ibis  thesis  addresses  the  problee  of  the  vertical  and 
horizontal  dispersion  of  the  forced  inertial  wave  energy 
fros  the  sized  layer.  The  forcing  readjusts  the  aean  flow 
as  veil  as  generating  anisotropic  inertial  notion.  The  cur¬ 
rent  aeter  observations  allow  detersination  of  the  scales 
and  the  energetics  of  the  force d,  anisotropic  wave  sotion 
through  use  of  spectrua  analysis  and  coaplex  demodulation. 
The  aodal  structure  is  exaained  by  fornulating  the  stura- 
Liouville  problea  and  then  perforaing  a  least  squares  fit 
between  the  eigenfunctions  and  the  deaodulated  time  series 


of  the  current  observations. 
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II.  HISTOBIC  hi  BEVIEW 


k  nonlinear  nodal  for  a  stationary  hurricane  (O'Brien 
and  Beid,  1967)  was  developed  to  sinulate  the  observations 
eade  in  the  wake  of  hurricane  Hilda  (Leipper,  1967).  The 
results  indicated  that  upvelling  regions  were  restricted  to 
the  area  enclosed  by  the  caxiiUB  wind  regise,  whereas  down- 
welling  occurred  outside  this  regiae  as  warn  water  was 
advected  away  froe  the  stcra  center. 

Cne  of  the  aost  proncunced  responses  of  the  ocean  to 
wind  forcing  is  the  generation  of  inertial  oscillations  in 
the  surface  aixed  layer.  Pollard  (1970)  siaulated  the  gen¬ 
eration  of  these  waves  using  a  t wo-diaensional,  wind- driven 
aodel  developed  by  Pollard  and  Hillard  (1970).  The  pre¬ 
dicted  aaplitudes  and  decay  rates  cf  the  forced  oscillations 
agreed  well  with  the  Woods  Hole  Oceanographic  Institution 
(W HOI)  'site  D'  observations.  Using  the  sane  data  set,  Pol¬ 
lard  (1980)  showed  that  under  relatively  strong  wind  condi¬ 
tions,  67  X  of  the  horizontal  kinetic  energy  (HKE)  in  the 
aixed  layer  was  found  near  the  inertial  frequency.  The  sub¬ 
sequent  downward  propagation  of  energy  was  estiaated  to  be 


of  the  order  of  10-*  ca/s  (  1  m/d) .  This  vertical  group 
velocity  was  too  saall  to  account  for  the  radiation  of  iner¬ 
tial  waves  and  the  subsequent  loss  of  energy  froa  the  sized 
layer.  The  corresponding  scales  of  the  inertial  wave  motion 
were  about  100  to  240  ■  in  the  vertical  and  hundreds  of 
kilometers  in  the  horizontal  (Table  I)  . 

TABLE  I 

Coapariscn  of  Inertial- Internal  lave  Paraneters 


Length 

Scales 

Horizontal  (km) 
Vertical  (a) 

Brooks 

Price 

Pol lard 

(1983) 

(1983) 

(1980) 

370 

480 

700-1700 

1000 

1000 

100-240 

Group 

Velocity 

Horizontal  (ka/d) 

23 

86 

0. 8-17. 

vertical  (a/d) 

60 

160 

0.03-3. 

|  Proa  Brccks  (1983) 

Osing  a  linear,  two- layer  acdel,  Seisler  (1970)  siau- 
lated  inertic-gravity  waves  as  part  of  the  baroclinic 
response  in  the  wake  of  a  translating  stora.  The  critical 
factcrs  governing  the  generation  of  these  waves  were: 

•  the  translational  speed  of  the  hurricane  nust  be 
greater  than  the  internal  wave  phase  speed;  and 

•  the  horizontal  scales  of  these  waves  sust  be  cosparable 
to  the  oceanic  Bcssby  radius  of  deformation. 


He  further  noted  that  the  ocean  response  to  a  aoving  hurri¬ 
cane  is  strongly  baroclinic. 

The  ocean  theraal  response  to  hurricane  forcing,  was 
investigated  nuaericaUy  by  iaposing  Ekaan  layer  dynaaics  in 
a  aired  layer  aodel  (Elsberry  et  a^.,  1976).  The  input  of 
energy  frca  the  wind  stress  was  forced  to  generate  inertial 
oscillations  and  cause  turbulent  airing  via  entrainaent. 
Op  welling  also  enhances  the  turbulent  airing  process  in  the 
upper  layers  of  the  ocean.  The  sain  result  was  that  advec- 
tion  doainated  the  theraal  response  near  the  storm  track  as 
opposed  to  the  heat  loss  froa  the  ocean  surface  to  the 
store. 

Strong  ataospheric  forcing  enhances  the  turbulent  airing 
process  in  the  surface  aired  layer  and  causes  it  to  deepen. 
As  the  aired  layer  continues  to  deepen,  the  theraocline 
begins  to  erode  as  water  froa  the  theraocline  is  entrained 
into  the  aired  layer.  However,  stratification  tends  to  sup¬ 
press  turbulent  airing  and  aids  in  the  creation  of  saall- 
scale  internal  waves  either  through  entrainaent  or 
Kelvin-Helaholtz  instabilites  (KB)  at  the  base  of  the  aired 
layer  (Pollard  et  §1. ,  1973  ;  Garwood,  1977). 


Below  the  nixed  layer,  ocean  current  variability  is 
linked  to  vertically  propagating  wave  groups  of  the  internal 
wave  field  (Kase  and  Olbers,  19  79)  .  The  actual  mechanisms 
for  the  vertical  transport  of  energy  from  the  wind  forced 
aixed  layer  tc  the  theraccline  for  the  generation  of  large 
scale  inertial-internal  waves  are  not  well  understood.  Per¬ 
haps  the  aost  effective  way  of  creating  these  large  scale 
inertial-internal  waves  is  through  Ekaan  suction.  Krauss 
(1972  a,b,  1976)  showed  that  a  horizontally  varying  wind 
stress  causes  a  mass  transport  9  0  deg.  to  the  right  of  the 
stress  in  the  Ekaan  surface  layer.  This  surface  layer  div¬ 
ergence  is  acccapanied  by  an  upward  displaceaent  of  the  iso- 
pycnals  or  upwelling  of  cccler  water  froa  below  which  tends 
to  suppress  the  turbulent  aixing  process.  During  the  relax¬ 
ation  of  the  wind,  these  isopycnals  are  displaced  downward  , 
which  contributes  to  the  creation  of  large  scale  inertial- 
internal  waves  at  the  base  of  the  aixed  layer. 

Curing  the  passage  of  hurricane  Belle  (Bayer  et  al . , 
1981),  ocean  current  and  temperature  measurements  were 
acquired  on  the  continental  shelf  of  the  Biddle  Atlantic 
Bight.  Their  analyses  indicated  that  aost  of  the  inertial- 
internal  wave  energy  was  ccntained  in  the  first  node  at  the 
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deeper  sites  (eater  depth  of  70  ■)  and  in  a  heavily  damped 


second  aode  at  the  shallow  sites  (eater  depth  of  50  m)  .  The 
variability  in  the  aodes  eas  not  attributed  to  the  spatial 
variability  of  the  wind  stress.  The  bottoa  slope  and  mean 
velocity  fields  significantly  altered  the  oceanic  response 
to  hurricane  Eelle. 

A  mixed  layer  model  (Garwood*  1977)  was  embedded  into  a 
multi-layer*  primitive  equation*  ocean  circulation  model 
(GCH)  by  Adamec  et  aj.*  (1980).  Using  this  GCH  model,  Hop¬ 
kins  (1982)  simulated  the  baroclinic  response  to  a  forcing 
pattern  similar  to  hurricane  Frederic  and  compared  the 
results  to  the  data  collected  in  the  wake  of  Frederic  (Shay 
and  Taaul*  1980).  The  model  predicted  the  inertial  response 
in  the  mixed  layer  quite  well.  However*  the  predicted  ocean 
current  response  in  the  subsurface  layers  was  less  energetic 
than  the  observations,  and  decayed  much  too  fast.  These 
discrepancies  between  the  model  simulations  and  observations 
are  due  in  part  to  the  imposition  of  the  rigid  lid  condition 
at  the  surface*  which  eliiinates  the  barotropic  aode.  Fur¬ 
thermore*  these  model  simulations  did  not  include  topogra¬ 
phy*  even  though  the  observations  are  from  a  region  where 
the  bottoa  topography  is  rugged  and  the  ocean  depth  is  less 
than  1  km  deep. 
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Recently#  ocean  measurements  were  aade  during  the 
passage  of  hurricane  Allen  in  the  western  Gulf  of  Mexico 
(Brocks#  1983) .  The  spatial  scales  of  the  inertial  wave 
were  different  froa  the  'site  D'  observations#  as  shown  in 
Table  I.  The  vertical  scales  in  the  ocean  are  generally 
auch  greater  under  hurricane  forcing  than  during  the  passage 
of  a  cold  front#  while  the  converse  is  true  for  the  horizon¬ 
tal  scales  cf  notion.  The  vertical  group  velocity,  which 
depends  on  both  the  wave  and  the  Brunt-Vaisala  frequencies, 
for  the  Allen  observations  exceeded  the  group  velocity  froa 
the  'site  O'  observations.  This  estinate  froa  the  Allen 
observations  is  nisleading  due  to  the  lack  of  upper  thermo- 
cline  and  aixed  layer  data.  Moreover#  the  wind  speed  is 
auch  greater  during  the  passage  of  a  hurricane  than  a  fron¬ 
tal  passage.  Therefore#  the  ancunt  of  turbulent  mixing  in 
the  aixed  layer  and  inertial  wave  excitation  should  be  auch 
aore  rapid  during  a  hurricane.  Other  features  measured  in 
the  Allen  response  included  topographical  dependence#  and 
the  vertical  phase  locking  during  the  first  few  inertial 
period  (IP)  following  the  stora. 

The  baroclinic  response  of  the  ocean  to  a  hurricane  was 
aodeled  by  (Erice#  1983)  using  a  multi-level#  inviscid 


model.  The  scales  and  energetics  of  the  simulated  inertial 
response  were  similar  to  mixed  layer  data  collected  by  a 


90 A&  data  huoy  during  the  passage  of  hurricane  Eloise.  The 
horizontal  and  vertical  scales,  as  given  in  Table  I,  wer9 
quite  large  in  comparison  to  the  thickness  of  the  thermo** 
cline  (200  a)  .  The  maxima  energy  of  the  inertial-internal 
wave  motion  mas  predicted  to  occur  at  a  distance  of  twice 
the  radius  of  maximum  winds,  hereafter  referred  to  as  the 
maximum  wind  regime,  which  was  roughly  80  km  for  Eloise. 
Further,  the  rate  of  vertical  energy  propagation  was  large 
in  comparison  to  both  the  illen  and  'site  D'  observations, 
and  accounted  fcr  the  depletion  of  energy  from  the  mixed 
layer.  Other  mechanisms,  such  as  KH  instability  (Pollard  et 
ai.,  1973  ;  Garwood,  1977)  and  turbulence  (Bell,  1978),  act 
to  remove  energy  from  the  mixed  layer.  However,  these  phe¬ 
nomena  were  neglected  as  mixing  was  not  explicitly  included 
in  the  model  formulation. 

Greatbatch  (1983)  modeled  the  nonlinear  response  of  the 
ocean  to  a  moving  storm.  The  major  results  were  that  the 
transition  between  upwelling  and  downwelling  zones  in  the 
oscillatory  wake  was  rapid,  and  that  nonlinearities  account 
for  the  displacement  of  the  maximum  response  to  the  right  of 
the  storm  track. 
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V-viv: 


III.  DATA  AN£  ANALYTICAL  METHODS 
A.  CURRENT  METER  DATA 

Ten  Aanderaa  RCM-5  current  meters  were  deployed  on  three 
moored  taut-wire  arrays  ic  depths  ranging  from  100  to  470  m 
in  the  northern  Gulf  of  Mexico.  These  current  meters  sam¬ 
pled  ocean  current  speed,  direction  and  temperatures  at  10 
minute  intervals.  Two  current  meter  arrays  (CMA2,  and  3) 
were  deployed  on  adjacent  sides  of  the  Desoto  Canyon  where 
bathymetric  ccntcurs  converge  to  form  the  head  of  the  canyon 
(see  Eig.  1)  .  The  other  scoring  was  deployed  closer  to  the 
coast  in  a  depth  of  about  100  a  of  water  where  the  isobaths 
are  nearly  parallel  to  tke  coast.  A  synopsis  is  given  in 
Table  II  cf  tke  storm  period  observations,  which  extends 
from  a  few  days  prior  to  hurricane  passage  (12  September)  to 
the  end  of  the  deployment  period  (mid-October)  . 

The  quality  cf  the  data  is  generally  good;  however,  the 
lengths  of  the  time  series  are  not  all  the  same.  For 
instance,  the  Savonius  rctors  were  eventually  lost  from  all 
current  meters  in  the  mixed  layer  due  to  the  large  current 
speeds.  These  large  nixed  layer  current  speeds  are  not 


corrected  for  rotor  pumping,  which  is  a  function  of  the 
mooring  design  as  well  as  the  large  direction  vane  of  the 
Aanderaa  current  meters  (Fofonoff  and  Ercan,  1967).  in 
addition  to  the  rotor  problems,  temperature  records  from  all 
mixed  layer  current  meters  were  unrecoverable  because  the 
ocean  temperature  exceeded  the  threshold  temperature  of  the 
thermistors  of  21.5  °C  .  Sea  surface  temperatures  325  km 
south-southeast  of  the  mooring  sites,  as  measured  by  NCAA 
Buoy  42003,  reached  28.8<>C  (Johnson  and  Renwick,  1981)  . 

B.  STORM  TRACK 

The  storm  track  is  based  on  the  best  position  data  on 
hurricane  Frederic’s  movement  through  the  Gulf  of  Mexico 
(Hebert ,  1979)  .  Because  of  the  intensity  of  Frederic,  recon¬ 
naissance  aircraft  constantly  monitored  the  storm.  Frederic 
increased  to  maximum  strength,  maximum  winds  or  minimum 
barometric  pressure,  80  to  130  km  west  of  the  cma  sites 
about  2100  GMT  12  September  (see  Fig.  1).  A  visible  photo¬ 
graph  of  Frederic  from  a  GOES  satellite  at  2001  gmt  12 
September  clearly  delineates  a  well  developed  eye  of  about 
40  to  50  km  in  diameter  (Figure  2) .  The  translational  speed 
of  the  hurricane  at  this  time  was  about  7  to  7.5  m/s  as  it 
approached  the  Gulf  Coast,  which  was  much  larger  than  the 


internal  va ve  phase  speed.  Thus,  an  inertio-gravity  wavs 
response  is  expected  in  the  wake  of  hurricane  Frederic 
(Geisler,  1970). 
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TABLE  II 

A  Synopsis  of  Stern  Period  Observations 
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C.  SIND  FIELD 

Bind  field  data  were  obtained  froa  shore  stations  at 
Mobile,  Alabaaa  and  Pensacola,  Florida.  These  records  indi¬ 
cated  that  the  surface  wind  speed  did  not  exceed  40  a/s. 
However,  ether  reports  suggested  that  wind  speeds  ranged 
between  48  to  58  a/s  as  Frederic  aade  landfall  (Hebert, 
1979).  These  discrepancies  are  attributed  no  local  boundary 
effects  and  the  non-representativeness  of  wind  data  col¬ 
lected  in  the  coastal  region.  Marine  winds  were  also  aeas- 
ured  by  NO  AA  Euoy  42003  (Jchnson  and  Renwick,  198  1).  Bind 
speeds  at  this  location  never  exceeded  35  a/s.  The  eye 
clearly  passed  over  the  buoy  as  indicated  by  the  ainiaua  in 
the  wind  field  as  the  direction  changed  froa  40  to  200  ° 
True  (Fig.  3),  with  a  corresponding  decrease  in  pressure  to 
959  aillibars  (ah). 

D.  VERTICAL  TEMPERATURE  GRADIENTS 

The  AXBT  data  collected  by  the  National  Hurricane 
Research  Division,  as  reported  by  Black  (1983),  are  used  for 
the  ccaputaticn  of  the  Brunt-Vai sala  frequency.  These  data 
were  collected  after  the  passage  of  Frederic  in  the  area  of 
current  aeter  array  deployaents  near  the  DeSoto  Canyon  (see 
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wind  speed  tie!  series  where  the  abscissa  depicts 
tiae  in  Julian  Oars  starting  an  1200  GHT  11  Sept, 
to  1500  GHT  13  sept.  1979.  tha  lower  panel 
represents  the  concurrent  wind  direction  tine  , 
series.  The  hatched  area  depicts  the  period  when 
the  data  buoy  was  in  the  eye  of  the  hurricane 
(Johnson  and  Renwiclc,  1981). 


The  vertical  density  profiles  are  computed  at  10  o 
intervals  using  the  AXBT  data  and  climatological  data  from 
N&VOCEANO  through  the  equation  of  state 

P  *  Pp(l-a(T  -  T0))  9  (1) 


where  p  is  the  density,  T  is  the  observed  temperature,  T« 
and  p.  are  the  reference  temperature  and  density  respec¬ 
tively  from  climatology,  and  a  is  the  thermal  expansion 
coefficient,  which  is  taken  to  be  0.0002/°C  .  some  error  is 
expected  due  tc  the  neglect  of  the  salinity  term.  However, 
an  examination  of  the  climatological  T-s  diagram,  for  the 
DeSoto  Canyon  area  indicates  that  density  variations  are  due 
to  temperature  rather  that  salinity  effects.  The  Brunt-Vai- 
sala,  H*  ,  frequency  is  computed  using  a  centered  finite 
difference  free  the  expression 
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where  g  is  the  acceleraticr  due  to  gravity  anddp/Az  is  the 
vertical  density  gradient.  Since  the  AXBT  data  only 
extended  to  about  250  a,  it  was  necessary  to  extrapolate  the 
vertical  teapexature  gradient  to  470  m.  The  temperature  is 
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The  lccal  variations  in  the  bottom  topography  at  the 
OeSotc  Canyon  (see  Pig.  1)  are  appproxioately  parallel  and 


normal  to  the  coast  at  CM12  and  CHA3.  Por  example,  north- 
south  flow  at  CHA2  is  in  the  crcss-shelf  direction,  but  it 


is  in  the  alcng-shelf  direction  at  CMA3.  Henceforth,  a 
Cartesian  coordinate  system  is  used  to  represent  the  hori¬ 
zontal  velocity  components  in  the  x  and  y  directions,  where 
eoticn  to  the  east  and  the  north  corresponds  to  +u  and  +  v 
current  velocities,  respectively. 

The  fundamental  time  scale  used  in  this  thesis  is  an 
inertial  pericd,  which  is  equal  to  24.  10,  24.25  and  24.47  h 
at  CHA  1,2  and  3  respectively.  In  the  following  calcula¬ 
tions,  the  Inertial  Period  (IP)  will  be  set  equal  to  24.4  h. 
The  intrinsic  length  scales  are  either  the  radius  of  maximum 
winds,  approximately  40  km  for  Prederic,  or  the  Rossby  radii 
of  deformation,  which  are  given  in  Table  III  for  the  baro- 
tropic  and  first  three  baxcclinic  modes.  These  deformation 
radii  are  based  cn  the  inertio-g ravity  wave  phase  speed  com¬ 
puted  from  the  Sturm-Liouville  problem,  which  will  be  dis¬ 
cussed  in  a  later  section.  For  convenience,  the  scales  in 


the  x  and  y  directions  will  be  referenced  as  cross  and 
along-track  scales  respectively,  because  the  storm  track  was 
nearly  perpendicular  to  the  coastline. 


TABLE  III 

Bossby  Radii  of  Deformation 


node 

No. 


Phase  Deformation 

Speed  Radii 

(m/s)  (km) 


68.0 

4.0 

2.2 

1.3 


950. 

56. 

30. 

18. 


F.  COMPUTATION  OF  SPECTRA 

The  energy  and  cross  spectra  are  computed  in  the  follow- 
ing  manner:  (1)  the  mean  is  computed  and  is  extracted  from 
the  time  series;  (2)  a  Tukey  data  window  is  applied  to  the 
data;  (3)  the  data  are  transformed  into  frequency  space  via 
the  Fast  Fourier  Transform  (FFT)  ;  and  (4)  the  transformed 
data  are  spectrally  averaged  over  bandwidths  which  vary  in 
width  with  frequency,  and  are  closely  tied  to  the  computa¬ 
tion  of  confidence  levels.  The  number  of  data  points  trans¬ 
formed  by  the  FFT  used  in  the  spectrum  analyses  does  not 
have  to  be  proportional  to  2n  ,  where  n  is  an  integer.  The 
general  FFT  is  advantageous  in  the  analysis  of  events 
because  of  their  transient  nature.  However,  in  the  analysis 
of  the  spectra,  equal  numbers  of  data  points  are  used  to 
avoid  leakage  and  smearing  of  energy  from  adjacent  frequency 


bands  (Otnes  and  Enochson,  1978)  .  Pinally,  the  higher  fre¬ 
quency  energy  abcve  the  Nyguist  frequency,  1/(2  A  t)  where 
A  t  is  the  sampling  interval,  is  elisinated  to  prevent  ali¬ 
asing  of  the  spectra. 

The  HKE  is  decoeposed  into  the  clockwise  (CW)  and  coun¬ 
terclockwise  (CC8)  rotating  components  using  rhe  rotary 
spectrum  method  outlined  in  Sonella  (1972)  and  nooers 
(1972).  Th9  main  results  of  this  analysis  are  the  rotary 
spectrum,  ellipse  stability  and  orientation,  and  rotary 
coefficient.  The  orientation  of  the  ellipse  indicates  the 
direction  of  the  horizontal  phase  velocity,  while  the  sta¬ 
bility  indicates  the  variability  in  the  ellipse  orientation, 
l  high  stability  index  indicates  that  the  wave  motion  is 
anisotropic,  cr  that  the  horizontal  direction  of  phase  pro¬ 
pagation  is  unidirectional.  Conversely,  a  low  index  of  sta¬ 
bility  implies  that  the  direction  of  phase  propagation  is 
random  or  isotropic.  The  rotary  coefficient  indicates  the 
amount  of  confidence  that  can  be  placed  in  the  rotary  spec¬ 
tral  estimate  and  th9  rotational  component  that  dominates 
the  rotary  spectrum.  For  instance,  a  -*>1.00  implies  that  the 
motion  is  polarized  in  the  CW  direction.  k  more  rigorous 
mathematical  treatment  is  given  by  Gonella  (1972)  and  Hooers 


G.  COMPLEX  OE EODULATION 


A  useful  method  for  determining  inertial  oscillations  in 
a  tiae  series  of  ocean  current  aeasurements  is  through  the 
use  cf  the  coupler  deaodulation  aethod  (Perkins,  1970).  The 
matheaatical  details  are  given  in  Appendix  A.  This  method 
can  be  inpleaented  either  by  linear  filtering  or  by  perform¬ 
ing  a  least  sguares  harmonic  analysis  on  the  time  series. 
The  linear  filtering  aethod  is  preferred  because  data 
points  are  lost  whenever  harmonic  analyses  are  used  to  form 
the  time-varying  amplitudes  and  phases  of  frequency  depen¬ 
dent  motions.  This  is  especially  useful  in  the  analysis  of 
events  when  all  data  points  are  needed  to  resolve  the  tran¬ 
sient  motions.  A  major  disadvantage  of  the  linear  filtering 
methcd  is  that  it  only  applies  to  discrete  frequencies, 
whereas  the  least  squares  method  can  treat  all  frequencies 
within  the  constraints  of  the  time  series. 

Initially,  the  data  are  band-pass  filtered  between  18 
and  30  h  to  limit  the  frequency  content  in  the  time  series 
to  near-inertial  motions.  These  filtered  data  are  multi¬ 
plied  by  the  trigonoaetric  arguments  of  the  inertial  fre¬ 
quency.  At  this  point,  some  higher  frequency  notions  have 
been  introduced  into  the  tiae  series  (Appendix  A)  and  these 


are  removed  by  low-pass  filtering  the  data  at  22  h.  The 


resultant  time  series  is  manipulated  to  form  the  instantane¬ 
ous  amplitudes  and  phases  cf  the  inertial  motion.  The  final 
set  of  time-varying  amplitudes  and  phases  provide  a  local 
harmonic  analysis,  or  nore  explicitly,  a  time  series  of 
spectral  estimates  for  inertial  motion.  The  filter  used  in 
the  demodulation  was  a  Lanzcos  square  taper  window  with  241 
weights  (Fig.  5)  . 

The  band-pass  filtered,  mixed  layer  velocity  components 
at  CHA3  are  excellent  examples  for  using  complex  demodula¬ 
tion  (Fig.  6).  As  the  current  velocities  increase  in 
response  to  the  storm,  the  amplitude  of  the  forced  inertial 
wave  increases  as  well.  The  inertial  wave  amplitudes  modu¬ 
late  the  increases  and  decreases  in  the  mixed  layer  veloci¬ 
ties.  After  the  maximum  velocities  occur,  roughly  2  IP 
following  hurricane  passage,  the  strength  of  the  currents 
decreases  linearly  with  e-folding  scales  of  about  4  ip. 
This  type  of  e-folding  behavior  is  very  common  in  amplitude- 
modulated  signals. 
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In  the  following  analyses,  the  data  are  divided  into  two 
periods:  1)  the  storm  period  observations  as  previously 
described  in  Table  II;  and  2)  the  quiescent  period  which 
encoapasses  the  period  frcn  the  start  of  the  data  to  a  few 
days  prior  to  the  stora.  This  chapter  focuses  on  the  iner¬ 
tial,  subinertial  and  superinertial  response  in  terms  of 
energetics  and  scales  of  action  to  hurricane  Frederic.  The 
quiescent  period  observations  and  results  are  given  in 
Appendix  B.  For  convenience,  the  terms  shallow  and  deep 
will  refer  tc  the  regicns  where  CMA1  and  CMA2,3  were 
deployed. 

A.  INERTIAL  BESECNSE 

Hurricane  Frederic  passed  abcut  80  Jem  west  of  the  DeSoto 
Canyon  at  2100  GMT  12  September  (JD  255)  1979.  The  current 
speeds  began  to  increase  at  070  0  GMT  or  approximately  14  h 
prior  to  the  tiae  of  closest  approach.  The  currents  in  the 
mixed  layer  at  CBA2  increased  to  a  maximum  within  an  IP  fol¬ 
lowing  the  passage  of  Frederic  (Fig.  7)  .  The  storm  caused 
the  horizontal  velocity  to  oscillate  with  periods  of  about 


24.3  h  and  to  rotate  in  a  CH  direction,  which  indicates  the 


Component  U  Component 


presence  of  inertial  waves  in  the  sized  layer.  These  waves 
then  decreased  in  amplitude  over  e-folding  scales  of  4  IP. 
The  storm  alsc  caused  the  u-component  to  flow  in  the  oppo¬ 
site  direction  from  that  cf  the  pre-storm  flow,  whereas  the 
v-component  increased  markedly  in  amplitude.  The  near-bet- 
tom  current  speeds  (shown  later)  began  to  increase  within  3 
h  after  the  sized  layer  response. 


Pigure  7. 


In  the  nixed  layer  at  CMA3 ,  the  storm  amplified  the 
pre-storm  currents  to  a  maximum  cf  135  cm/s  (Fig.  8).  As  at 
CM A2 ,  this  increase  in  the  currents  is  a  manifestation  of 
inertial  wave  excitation.  However,  the  inertial  response  at 
CMA3  was  different  from  CMA2  in  that  the  kinetic  energy  of 
the  waves  was  larger  and  persisted  longer  (approximately  2 
IP)  .  The  energy  of  these  waves  then  decreased  over  e-fold¬ 
ing  scales  of  4  IP.  The  increase  in  near-bottom  current 
(shown  later)  occurred  within  4.5  h  following  the  mixed 
layer  response. 

These  horizontal  differences  of  the  inertial  response  in 
the  mixed  layer  between  CBA2  and  CMA3  are  observed  in  the 
amplitude  changes  of  the  HKE  (Fig.  9)  .  The  initial  increase 
in  kinetic  energy  at  the  two  arrays  was  similar.  The  maxi¬ 
mum  in  kinetic  energy  occurred  first  at  CMA2  and  led  the 
mixed  layer  maximum  at  CM  A3  by  about  3  h.  The  differences 
in  HKE  levels  were  quite  pronounced.  For  example,  the  mixed 
layer  HKE  level  was  roughly  60  %  cf  that  found  in  the  mixed 
layer  at  CBA3 ,  even  though  CM A2  was  closer  to  the  eye  of  the 
storm.  This  large  difference  over  a  distance  of  34  kn  may 
be  due  in  part  tc  the  stronger  mean  current  at  CM A2,  which 
nay  have  reduced  the  amplitude  of  the  inertial  waves.  The 
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Figure  8.  Sizilar  to  Fig.  7,  except  at  CMA3. 


direction  of  tie  nan  flow  and  the  horizontal  phase  propaga¬ 
tion  of  the  inertial  waxes  (shown  below)  roughly  coincided, 
which  suggests  that  there  say  have  been  some  interaction 
between  then.  Following  these  naxina,  the  energy  decayed 
with  an  e-folding  scale  cf  4  IP  at  both  arrays.  Forced 
inertial  waxes  in  the  sized  layer  that  xary  in  space  and 
tise  are  one  sechaniss  to  drixe  inertio-graxity  waxes  in  the 
subsurface  layers  (Krauss,  1972a, b). 


The  vertical  differences  in  the  inertial  response  were 
even  aore  dramatic.  In  the  deep  thermocline  (251  m)  at 
CHA3,  the  increase  in  the  currents  was  sensed  within  2  h 
follcving  the  sized  layer  response  (Fig.  10).  The  initial 
mazieuo  occurred  within  an  IP  (JD  256)  after  the  winds 
relax,  which  was  followed  by  inertial  oscillations  in  the 
horizontal  currents.  In  the  presence  of  stratification,  the 
particle  actions  are  elliptical  rather  than  circular  as 
observed  in  the  mixed  layer.  The  amplitude  of  inertio-grav- 
ity  waves  decreased  slightly  following  the  initial  response. 
About  10  IP  later,  the  current  increased  to  a  secondary  max* 
imua  that  nearly  exceeded  the  initial  maximum.  The  u-velo- 
city  essentially  defined  an  amplitude-modulation  envelope 
that  collapsed  and  expanded  in  time  and  persisted  for  21  IP 
following  the  initial  maxisum.  The  near-bottom  currents  at 
CHA3  (shown  later)  behaved  in  a  similar  manner  except  that 
the  secondary  maximum  occurs  6  IP  following  the  storm. 

The  vertical  changes  in  the  currents  at  CHA3  are  even 
aore  apparent  in  the  progressive  vector  diagrams  (PVD)  , 
which  are  Lagrargian  representations  of  the  Eulerian  meas¬ 
urements.  That  is,  the  continuous  current  measurements  are 
integrated  with  respect  to  time  to  approximate  the 


In  the  thermcdine  (Fig.  lib),  the  store  amplified  the 
freely  propagating  inertic-gravity  wares.  However,  the 
dimensions  of  the  ellipses  are  distorted  because  of  the 
increases  in  the  mean  flew  between  JO  257  and  261.  This 
period  was  immediately  followed  by  decreases  in  the  lean 
flow  until  JC  265,  when  the  inertio-gra vity  response  was 
more  pronounced.  The  period  of  increasing  and  then  decreas¬ 
ing  mean  flow  corresponds  to  the  approximate  modulation 
envelope  and  lag  period  of  10  IP  between  JD  256  and  266 
(Fig.  10).  after  JD  266,  the  modulation  envelope  started  to 
collapse  as  the  inertio- gravity  waves  decayed  in  amplitude 


with  an  e-f elding  scale  of  about  4  IP. 

The  vertical  propagaticn  of  the  energy,  which  is  a  con- 
seguence  of  the  generation  of  inertio-gravity  waves  in  the 
thermccline,  is  indicated  by  the  vertical  differences  in  the 
HKE  between  the  mixed  layer  and  the  thermocline  (Fig.  12). 
In  the  mixed  layer,  the  zaximua  HKE  occurred  approximately 
one  IP  following  the  storm  and  persisted  over  the  next  2  ip. 
Subseguently,  a  relative  maximum  was  reached  in  the  deep 
thermocline  atcut  2  h  after  the  surface  maximum,  and  it  was 
only  5  to  8  X  of  the  maximum  HKE  in  the  mixed  layer.  Based 
only  on  these  initial  phase  differences  between  the  velocity 


W 


components  at  the  two  levels,  the  vertical  propagation  of 
energy  would  he  about  3  ca/s  (2.6  km/d)  .  The  thermocline 
HKE  then  reaained  nearly  constant  over  the  next  4  to  6  IP. 
On  JD  262,  the  theraocline  HKE  gradually  began  to  increase 
as  the  vertically  propagating  energy  from  the  mixed  layer 
accuaulated.  Using  this  10  IP  lag  period,  the  vertical  pro¬ 
pagation  of  energy  is  about  0.03  ca/s  (26  m/d)  .  The  first 
estiaate  of  2.6  ka/d  is  about  1  tc  2  orders  of  magnitude 
greater  than  the  Eloise  and  Allen  observations  (see  Table 
I)  .  However,  the  second  estiaate  of  '26  m/d  is  much  more 
comparable  to  the  Allen  estimates.  The  propagation  of 
energy  will  be  addressed  below  using  the  dispersion  relation 
of  inertio-gravity  waves.  The  phase  propagation  was  clearly 
upward  over  the  first  10  IE  following  the  storm.  Over  the 

subsequent  10  IP,  the  propagation  was  not  as  clear  as  phases 
0 

in  the  nixed  layer  became  erratic,  presumably  due  ro  the 
decreasing  energy  levels  or  spectral  broadening. 

The  inertial  peak  dcainated  the  HKE  spectra  at  both 
locations  and  at  all  depths.  The  change  in  the  inertial  HKE 
between  the  surface  and  bcttoa  layers  was  about  an  order  of 
magnitude,  while  the  energies  between  mid-depth  and  interme¬ 
diate  levels  did  not  change  significantly  at  CHA3  (Fig.  13)  . 


Although  the  inertial  peak  was  prominent,  the  variability  in 
spectral  shapes  was  due  tc  the  differences  in  the  background 
continuum  of  internal  waves  which  varies  with  depth  (Fu, 
1960)  and  bandwidth  averaging. 

The  HKE  spectra  are  decomposed  into  CR  and  CCW  rotating 
components  to  clarify  the  rotational  characteristics  of  the 
inertial  wave  motion.  The  spectral  peak  near  the  inertial 
frequency  was  shifted  about  3-6 X  higher  than  the  inertial 
frequency  and  it  dominated  the  CR  spectra  calculated  from 
the  mixed  layer  velocity  measurements  at  CNA3  (Fig.  14) . 
The  rotary  spectrum  estimates  for  the  inertial  period 
motions  are  given  in  Table  IV.  The  CR,  inertial  HKE  vari¬ 
ability  estimates  are  generally  1  to  2  orders  of  magnitude 
larger  than  the  CCR  rotating  motion.  In  comparison,  the 
subinertial  (periods  of  3  IP)  and  superinertial  (periods  of 
12  h)  CR  spectral  estimates  are  more  than  an  order  of  magni¬ 
tude  less  than  the  inertial  motion  (Tables  VII  and  VIII)  . 

At  all  levels  at  the  deep  region,  the  positive  rotary 
coefficients  approached  unity,  which  indicates  that  the 
waves  were  polarized  in  the  CR  direction  (Table  IV) .  The 
direction  of  the  semi-major  axes  of  the  inertial  ellipses  at 
CBA2  rotated  CR  with  depth  except  near  the  bottom.  The 


Figure  13. 


Horaalized_HKE  Spectra  froa_CHA3  at  21  a 

a  (dashed-dotted) 
>resents  the  diurnal 


(solid),  251  a  (dashed),  437 

ind  457  a  (dotted).  K.  reprt  _  _ _ 

tide  frequency,  f  the  inertial  frequency  and  H. 
the  seai-diurral  tide  frequency. 
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Figure  14.  Hired  Layer  Serialized  Rotary  Spectra  from  CHA3 
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TABLE  IV 

Normalized  Rotary  Spectrua  Analysis  at  Inertial  Frequency 

Normalized 

Rotary 

Spectra  Ellipse 

Eepth  CU  CCU  Dir.  Rotary 

(a)  (cm/s)  2/cph  Stab.  (Dec.)  Coef f . 


a) 

CM! 

49 

6.7x10+3 

1.3x1  0+3 

0.81 

18 

♦  0.68 

64 

2.5x10+3 

1. 4x1 0+3 

0.20 

65 

♦  0.62 

b) 

£M2 

19 

1.8x10+3 

3.5x1  0+3 

0.90 

60 

♦  0.96 

179 

4.6X10+  + 

3.  7x1  0+2 

0.95 

80 

+  0.95 

324 

2.8x10+  + 

1.3x1 0+2 

0.30 

60 

♦  0.96 

s) 

CM3 

21 

4.0x10+3 

2. 0x1 0+3 

0.87 

50 

♦  0.98 

251 

1.4x10+3 

6. 1x1 0+2 

0.63 

90 

♦  0.98 

437 

1.2x10+3 

6.3x10+2 

0.50 

150 

♦  0.98 

457 

6.2x10+  + 

5.0x1 0+2 

0.77 

160 

♦  1.00 

direction  of  the  near-bottoa  ellipse  rotated  CCU  with  depth 
between  the  tteraccline  (179  a)  and  near  the  bottom  (324  m)  . 
The  stability  of  the  ellipses  indicated  that  the  inertial 
moticn  near  the  bottom  at  CHA2  was  isotropic.  In  contrast, 
inertial  oscillations  were  strongly  anisotropic  at  all  of 
the  ether  depths  in  the  deep  region.  The  inertial  wave 
response  at  CHA1  was  markedly  less  intense  than  at  CMA2  or 
CH A3 ,  bat  it  was  polarized  in  a  CU  sense.  The  reasons  for 
the  semi-diurnal  tide  at  CHA1  being  much  more  energetic  than 
the  inertial  motion  is  explained  below. 


The  vertical  and  horizontal  coherences  are  computed  to 
estimate  the  spatial  scales.  The  spatial  scale  L  is  com¬ 
puted  from  the  cross -spectral  phase  relationship  using 


t  360  .  . 
L  =  Air  Ad 


where  Ad  is  either  the  vertical  or  horizontal  distance 
between  the  current  meters  and  A4>  is  the  phase  difference 
between  similar  velocity  components  (Pollard,  1980)  .  This 
analysis  is  restricted  to  the  CHA2  and  3  due  to  the  absence 
of  an  energetic  inertial  response  at  CM&1  (Table  V) .  Esti¬ 
mates  of  vertical  scales  are  more  accurate  than  the  horizon¬ 
tal  scale  estimates  because  the  vertical  spacing  of  the 
current  meters  below  the  mixed  layer  is  not  the  same  at  the 
two  arrays.  These  scales  calculated  from  the  entire  storm 
time  series  will  be  compared  belcw  with  the  scales  calcu¬ 
lated  from  the  instantaneous  phases. 

The  vertical  cross- spectral  analyses  (Table  V)  indicate 
that  the  inertial  response  in  the  mixed  layer  led  the  near- 
bottom  response.  The  normalized  cross- spectrum  variance 
between  the  vertical  levels  was  about  a  factor  of  two 
greater  at  CHA3  than  at  CHI2.  The  vertical  scales  computed 
from  (3)  are  about  360  and  54  0  m  at  CHA2,3  respectively. 


That  is,  the  inertial  motion  had  scales  of  the  order  of  the 
water  depths  at  both  arrays  and  was  coherent  at  the  95  X 
confidence  level. 


TABLE  T 

Vertical  and  Horizontal  Coherencies  at  the  Inertial 

Frequen  cy 


Ar-ay  Variables 


Nortalized 
Cross-Spectra 
ZASiH SJZSil 


Squared  Phase 

Coherence  (deg. ) 


Vertical 


2 

u  (s) 

,u  <b] 

1  3.2x10++ 

0.62 

-50 

2 

V  (s) 

,v  <b' 

|  2. 7x10++ 

0.60+ 

-70 

3 

u  |sj 

,  u  (b 

6.5x10++ 

0.70 

-55 

3 

V(s) 

» v  (bj 

1  6. Cx10++ 

0.80 

-50 

Horizontal 

3,2 

u  <s] 

) 

1.9x10+8 

0.98 

-45 

3,2 

V  (s| 

1 

1.6x10+s 

1.00 

-55 

3,2 

u  (b 

) 

2.3x10++ 

0.98 

-20 

3  2 

v  (b| 

1 

2. 1x10++ 

0.95 

-30 

•*  net  significant  at  95 
*  spectra  assumes 
js|  =  surface 


X  confidence 


bottem 


mes  an  average  inertial  frequency  c 
layer  (CM A 2  *  19  a  ;  CMA3  =  21  mf 
layer  (CM A2  =  324  a  ;  CHA3  =  457  m) 


of  0.  0412  eph 


The  instantaneous  phase  differences  of  the  velocity  com¬ 
ponents  between  the  two  records  can  also  be  used  to  estimate 
the  vertical  scales  (Fig.  15).  Initially,  the  estimated 
vertical  scales  at  CHA3  were  much  greater  (typically  800  to 
900  a)  than  the  water  depth,  but  after  JO  257  the  estimates 
were  consistently  of  the  order  of  the  water  depth.  The 


large  vertical  scales  observed  initially  agree  with  these 
estimated  froa  the  Eloise  siaulations  by  Price  (1983)  and 
the  lllen  observations  by  Erooks  (1983).  At  CMA2 ,  the  ver¬ 
tical  scales  had  siailar  trends  and  converge  to  scales  of 
the  erder  of  the  water  depth. 

The  horizcntal  cress- spectrum  estimates  (Table  V)  indi¬ 
cate  both  velocity  coaponents  were  coherent  in  the  mixed 
layer.  The  corresponding  scales,  as  estimated  from  (3),  are 
20  and  35  km  in  the  east-west  and  north-south  directions 
using  an  array  separation  of  34  km.  Initially,  the  horizon¬ 
tal  scales  were  set  by  the  atmospheric  forcing,  but  these 
estimates  nay  be  somewhat  biased  because  the  entire  stern 
tine  series  was  used  in  the  computation  of  the  spectra.  The 
cross-spectrum  variance  between  the  two  near-bottom  records 
was  about  cne  order  of  magnitude  less  than  in  the  mixed 
layer.  However,  these  values  are  significant  at  the  95  % 
confidence  level  as  indicated  by  the  coherence  squared.  The 
respense  at  CBA2  leads  the  response  at  CHA3  by  about  2  h. 
The  horizontal  scales  near  the  bottom  calculated  froa  (3) 
are  20  and  35  km  in  the  east-west  and  north-south  direc¬ 
tions,  respectively.  The  differences  in  the  scales  are  due 
to  the  larger  scales  set  up  by  the  hurricane  in  the 


Figure  15.  Vertical  Scales  at  CIU3  computed  from  the 

instantaneous  phase  differences  between  the  0 
velocity  components  at  21  a  and  457  m. 


along-track  (ncrtb-south)  than  in  the  cross-track  (east- 
west)  directions.  These  estimates  of  scales  nay  be  somewhat 
misleading  because  of  the  instrument  depth  difference  (135 
m)  ,  although  they  agree  well  with  mixed  layer  estimates. 

The  instantaneous  phase  differences  between  the  mixed 
layer  currents  indicate  that  the  horizontal  length  scales 
were  set  bj  the  storm  initially  (Fig.  16).  The  scales 

decreased  rapidly  and  approached  constant  values  of  60  and 
100  km  in  the  east-west  and  north-south  directions,  respec¬ 
tively.  These  scales  in  the  mixed  layer  are  roughly  1-2 
times  the  baroclinic  Bossby  radius  of  deformation  for  the 
first  mode  (56  km).  The  estimates  of  the  horizontal  scales 
of  inertic-gravity  waves  in  the  tharmocline  are  computed 
from  the  tottcs  and  mid-depth  instantaneous  phase  differ¬ 
ences  in  the  0-velocity  component  at  CHA2  and  CHA3,  respec¬ 
tively  (Fig.  17) .  The  east-west  and  north-south  scales  in 
the  thermocline  immediately  after  hurricane  passage  are 
roughly  the  same  as  in  the  mixed  layer.  After  4  ip,  the 
scales  revert  tc  the  pre-storm  scales  of  25  and  50  km  in  the 
east-west  and  north-south  directions,  respectively.  seme 
caution  is  advised  in  accepting  these  estimates  due  to  the 
large  depth  difference  of  73  a.  In  a  later  section,  it  is 


suggested  that  the  energy  was  contained  primarily  in  the 
barotropic  and  first  two  baroclinic  modes.  Hence,  a  depth 
difference  of  13  a  may  be  acceptable  within  the  deeper  ther- 
aocline  because  of  the  relatively  large  displacement  of  the 
isotherms,  the  large  vertical  scales,  and  the  amount  of 
energy  contained  in  the  barotropic  mode. 

The  estimates  of  vertical  energy  propagation  speed  are 
computed  only  cn  the  basis  of  the  lag  times  between  the  sur¬ 
face  and  subsurface  response  (Table  VI).  In  the  first  case, 
the  response  is  felt  at  the  bottom  of  the  water  column 
within  3  to  4.5  h  after  increases  in  the  mixed  layer  are 
detected.  The  initial  response  corresponds  to  a  vertical 
scale  of  roughly  900  a  (Pig.  15).  This  lag  of  3  to  4.5  h 
results  in  a  vertical  energy  propagation  speed  of  2.6  ka/d 
which  is  roughly  1  to  2  orders  of  magnitude  greater  than  the 
values  calculated  by  Brooks  (198  3)  and  Price  (1983).  The  10 
IP  lag  in  the  second  maximum  observed  at  251  a  represents  a 
vertical  energy  propagation  speed  of  26  m/d  when  the  verti¬ 
cal  scale  is  about  500  a.  This  vertical  propagation  of 
energy  agrees  quite  well  with  the  estimates  derived  from  the 
Allen  observations.  However,  there  is  also  a  second  maximum 
observed  near  the  bottom  about  6  IP  following  the  storm. 


Consequently,  the  vertical  (Cj, )  and  horizontal  (CjW  ) 


group 

velocities  are 

calculated  froa 

the 

following  expres 

sions 

(Brooks,  1983) 

Cgz 

aN2  2 

-  i  ^2  tan  9 

7 

(4) 

CgH  = 

0  N2  2- 

k  ^2  tan  9 

7 

(5) 

where  <r  is  the  observed  frequency,  a  and  k  are  wavenuabers 
in  the  vertical  and  horizontal  cross- track  directions,  and 
tan0  *  k/a.  The  wavenumbers  a  and  k  are  equal  to  2  v /Laand 
2  it  /Irrespectively,  where  La  and  L„  represent  the  vertical 
and  horizontal  scales  as  estimated  previously.  The  assigned 
values  to  the  variables  in  (4)  and  (5)  are 

<r  *  7.54  x  10”s  rad/s, 

S2  *  3.  30  x  1  C-s  rad  */s 2, 
a  >  6.  98  x  10-2  rad/m, 
m  «  1.16  x  10-2  rad/m, 
k  a  1.05  x  10-*  rad /a. 

The  initial  vertical  group  velocity  is  1.25  ka/d  as  com¬ 
puted  froa  the  dispersion  relation  (4)  with  an  initial 


scale  of  900  a  (Table  VI) 


This  estimate  is  about  half  of 
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£ 

Lag  Time 
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.06 
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Group 

Velocity 

l 

P 

Dispersion 

(km/d) 

80.  30. 

23. 

86. 

*  From  Brook 

s  (1983) 

the  one  computed  based  on  the  initial  lag  tine  and  is  still 
an  order  of  magnitude  larger  than  the  estimates  computed 
from  the  Eloise  simulations.  After  a  few  IP,  the  vertical 
group  velocity  is  about  0.27  km/d,  which  is  of  the  order  of 
the  Eloise  estimates  of  about  0. 16  km/d  (Price,  1983)  .  Fur¬ 
thermore,  this  estimate  is  about  one  order  of  magnitude 
larger  than  the  value  of  0.02  km/d  that  Brooks  (1983)  esti¬ 
mated  from  tie  Allen  observations.  Clearly,  the  initial 
storm  forcing  causes  energy  to  propagate  vertically  at  a 
significantly  higher  rate  than  after  a  few  IP  following  the 
storm. 


v 


The  horizontal  propagation  of  energy  is  a  much  more 


rapid  than  the  vertical  propagation  (Table  71) .  The  group 
velocity  estimate  for  the  Frederic  observations  is  initially 
80  ki/d;  and  it  decreases  to  about  30  km/d  a  few  IP  follow¬ 
ing  the  storm.  This  difference  is  due  to  the  tan2  depen¬ 
dence  on  m  as  previously  defined.  The  initial  horizontal 
group  velocity  estimate  agrees  with  the  Price  (1983)  esti¬ 
mate  whereas  the  lower  value  agrees  with  the  Brooks  (1983) 
estimate.  The  horizontal  phase  speeds  of  the  inertio -grav¬ 
ity  waves  in  the  wake  of  Frederic  are  about  70  cm/s  with  a 
vertical  phase  speed  of  roughly  1  cm/s. 

In  summary,  inertial  waves  maintained  a  constant  ampli¬ 
tude  over  1  to  2  IP  following  the  passage  of  Frederic. 
These  oscillations  in  the  mixed  layer  then  relaxed  over  an 
e-folding  scale  of  about  4  IP.  In  the  subsurface  layers, 
the  forced  inertio-gravity  waves  decayed  over  similar  tine 
scales  after  the  appearance  of  a  secondary  maximum.  The  lag 
between  the  initial  impulse  of  energy  and  the  second  subsur¬ 
face  maximum  varied  over  depth.  For  example,  this  period 
was  about  10  IP  at  251  a,  but  it  was  only  6  IP  at  457  a. 
These  lag  values  define  a  modulation  envelope  which  evolves 
in  time.  If  vertical  propagation  of  energy  from  the  mixed 


layer  alone  accounts  for  these  secondary  maxima,  why  do  the 
maxima  occur  later  in  the  deep  theroocline  than  near  the 

i 

bottca  ?  ether  mechanisms  besides  dispersion  of  energy  from 
the  aixed  layer  aust  be  responsible.  For  example,  nonlinear 
resonant  interactions  among  the  inertio-gravity  waves  could 
be  responsible  (HcComas  and  Bretherton  ,  1977).  Analyses  of 
these  interactions  are  beyend  the  scope  of  the  present  work. 

The  initial  vertical  scales  of  the  inertio-gravity  waves 
in  the  wake  of  Frederic  are  approximately  equal  to  the 
Eloise  and  Allen  estimates.  Hcwever,  after  a  few  IP  the 
vertical  scale  decreased  to  about  500  m  or  approximately  the 
water  depth.  The  horizontal  scales  estimated  from  the  Fred¬ 
eric  observations  differ  greatly  from  the  Allen  data  and 
Eloise  simulations.  The  horizontal  scales  during  Allen  were 
four  times  the  mixed  layer  horizontal  scales  observed  immed¬ 
iately  after  the  passage  of  Frederic.  This  difference 
increased  to  a  factor  of  about  eight  approximately  4  ip 
after  Frederic.  These  discrepancies  may  be  due  to  the  rug¬ 
ged  fcettos  topography  in  the  DeSoto  Canyon  region. 


B.  SUBINEBTIAI  BESPOHSE 

The  forcing  by  hurricane  Frederic  also  altered  the  "mean 


flow"  and  the  lower  frequency  motions  (periods  greater  than 
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3  IP).  Tbs  variability  and  response  of  the  mean  flow  to  the 
forcing  is  depicted  in  progressive  vector  diagrams  (PVD). 
The  ocean  current  data  were  low-pass  filtered  with  a  half 
power  point  at  26  h  to  isolate  the  forced  response  on  the 
lower  freguency  variability. 

The  mean  flow  in  the  sized  layer  at  CM  A  2  was  directed 
towards  the  east  prior  to  the  stcrm  (Fig.  18a).  During  the 
passage  of  hurricane  Frederic  near  JD  256,  there  was  only  a 
small  deflection  in  the  mean  flow.  The  subsequent  horizon¬ 
tal  displacements  indicate  that  the  mean  flow  was  about  26 
kn/d,  or  almost  twice  the  mean  current  during  the  pre-stcrm 
period.  The  mean  flow  at  179  m  (not  shown)  at  CHA2  also 
increased  towards  the  east. 

The  bottcm  PVD  at  C3A2  indicates  that  the  mean  flow  at 
rhat  level  was  acre  or  less  constant  (Fig.  18b).  However, 
the  direction  of  this  flow  was  opposite  to  that  observed  in 
the  mixed  layer.  Consequently,  the  bottom  PVD  was  nearly  a 
airrcr  image  of  the  surface  PVD,  although  the  bottom  flow 
was  slower  than  that  of  the  mixed  layer.  There  was  also  a 
superposed  lower  frequency  motion  of  about  3  IP  near  the 
bottom,  particularly  during  the  period  between  the  initial 
and  secondary  maxima  of  HKE  (first  6  IP  in  Fig.  9).  The 
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lower  frequency  oscillations  may  be  due  to  coastally  trapped 
waves,  which  were  enhanced  by  the  storm,  since  waves  of  the 
same  period  were  observed  during  the  quiescent  period 
(Appendix  E) .  These  waves  can  be  altered  by  both  bottom 
topography  and  continuous  density  stratification  (Wang  and 
Mooers,  1976).  A  study  of  the  mechanisms  by  which  the  mean 
flow  is  altered  by  these  waves  and  the  rugged  terrain  is 
also  beyond  tbe  scope  of  this  thesis. 

At  CHA3,  the  surface  mean  flow  was  southward  pricr  to 
the  storm  (Fig.  19a).  At  the  tine  of  storm  passage,  near  JD 
256,  the  surface  mean  flow  was  deflected  outward  from  the 
storm  center.  Between  JD  257  to  269,  the  mean  flow  acceler¬ 
ated  to  about  13  km/d  and  then  decelerated  to  about  6  km/d 
towards  the  south.  Prior  to  the  storm,  the  mean  flow  near 
the  bottom  at  CMA3  was  northward  following  the  local  bottom 
topography  (Fig.  19b)  .  There  was  an  eastward  (up  the  canyon 
slope)  deflection  late  on  JD  256.  After  completing  a  loop, 
the  mean  flow  veered  towards  the  northwest  (down  the  canyon 
slope).  Between  JD  254  to  263,  the  mean  flow  near  the  bct- 
tom  was  considerably  faster  than  either  earlier  or  later 
times.  Again,  the  mean  flow  near  the  bottom  was  reversed 
relative  tc  that  in  the  mixed  layer. 


Figure  18. 


Lew-pass  Filtered  PV C  at  CMA2.  The  *  represen 
a  tiae  interval  of  48  h  starting  on  JD  248  for 
a)  mixed  layer  (19  a)  and  b)  bofton  Layer  (324 
■  ) . 
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Figure  19.  Lev-pass  Filtered  PVE  at  Ci!A3.  The  *  represents 
a  ti?e  interval  of  48  h  starting  on  JD  245  for 
at  mixed  layer  (21  a)  ana  b)  bottom  Layer  (45v 

■  )  • 


14.5°c  in  the  theraocline.  At  CHA3r  the  temperature  at  251 
m  started  to  increase  about  0600  GHT  on  JD  255  and  continued 
to  increase  over  the  next  15  h  due  to  the  strong  northward 
current.  The  temperatures  then  decreased  by  2°c  due  to  a 
strong  southward  current.  Inertial  fluctuations  of  abcut 

After  this  time. 


3<>C  persisted  ever  the  subsequent  8  IP 


the  variations  underwent  a  slight  increase  in  frequency 
which  corresponded  to  about  the  tine  of  the  secondary  maxi¬ 
mum  in  the  horizontal  velocity  components.  For  the  remain¬ 
der  of  the  record,  the  temperature  variations  were 
superposed  on  a  stationary  mean. 

Ecttoa  temperature  variability  at  CSA2  is  a  clear  exam¬ 
ple  cf  advective  processes  that  are  associated  with  lower 
frequency  moticns  (Fig.  21).  On  JD  254  (11  Sept  ember)  ,  a 
temperature  decrease  of  1.3<>C  over  a  16  h  period  marked  the 
beginning  cf  the  forcing  period  as  cool  water  was  advected 
from  the  DeSoto  Canyon.  The  temperature  then  increased  by 
3.70  c  over  the  next  25  h  as  warm  water  was  advected  by  the 
southward  velocity  component  which  was  downslope  in  this 
part  of  the  canyon.  As  the  winds  relaxed,  cooler  water  was 
again  advected  into  the  area.  Around  JD  260,  the  tempera¬ 
ture  dropped  ty  4oc  ever  an  11  hour  period.  This  drastic 
change  in  temperature  was  due  to  the  advection  of  cold  water 
by  these  subinertial  waves  of  3  IP  period.  From  the  orien¬ 
tation  of  the  ellipses  (Table  ?II)  ,  the  direction  of  the 
wave  propagation  was  northward  (upslope,  away  from  the 
Desoto  Canyon).  These  waves  set  up  cross-shelf  oscillations 
(relative  to  CHA2)  in  the  north-south  direction.  After 
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these  large  variations,  near-inertial  oscillations  in  the 
teaperature  were  superposed  on  a  gradual  warning  trend. 

TABLE  VII 

Normalized  Rotary  Spectrun  Analysis  at  Subinertial  Frequency 

Moraalized 

Rotary 


Cepth 

Spectra 

Cl  CCl 

Elli  pse 

Dir. 

Rotary 

<B) 

(sa  zs) 

Stab. 

(Deq.) 

Coeff . 

324 

4.5x1 0+ 3  3.0X10+3 

0.65 

10 

♦  0.20 

457 

4.7x10+3  1.0x10+3 

0.50 

10 

♦  0.98 

The 

potential  energy  spectra 

were  computed 

using  the 

teaperature  tine  series  and  the  vertical  temperature  gradi¬ 
ents  froa  the  AXET  data  collected  by  Black  (1983)  (Appendix 
A)  .  The  nest  prcainent  peak  in  the  potential  energy  spectra 
at  CHA3  is  near  the  inertial  frequency  (Fig.  22).  The  most 
energetic  peaks  at  the  seai-diurnal  and  subinertial  frequen¬ 
cies  are  near  the  bottoa.  potential  energy  variability  in 
the  bottoa  layers  at  CHA2  (not  shown)  is  significant  at  the 
95  %  confidence  levels  only  for  periods  of  about  3  IP. 

Froa  rotary  spectrun  analyses  at  periods  of  3  IP  near 
the  bottoa:  (1)  the  Cl  and  CC1  rotating  components  at  CHA2 
were  nearly  equal,  with  just  a  slight  tendency  towards  Cl 
polarization;  and  (2)  the  notion  at  CHA3  was  polarized  in  a 
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Pigure  22. 


pcteptial  ?nergy  Spectra,  froa  CSA3 
lace  a 


isctfcera  <3 
(sclid)  ,  4 


._.jents  observed  at 
(dashed)  and  457 


|E  th. 


and  4&7  a  (dotted)  .  K. 
.•presents  the  diurnal' tide  frequency,  f  the 
Lnertial  frequency  and  Hx  the  seai-diumal  tide 
frequency. 
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CW  sense  (Table  VII)  .  Furthermore,  the  phase  propagation 
for  these  longer  period  waves  was  northward,  and  anisotropic 
as  indicated  ky  the  ellipse  orientation  and  stability. 

The  cross-spectrum  between  the  temperature  and  north- 
south  velocity  records  at  the  bottoe  was  also  dosinated  by 
the  peak  at  about  3  IF  (Fig.  23).  The  corresponding  phase 
and  coherence  estimates  are  100  deg.  and  0.98  respectively, 
with  the  velocity  leading  the  teaperature  signal.  Hence, 
these  lover  frequency  oscillations  in  the  cross-shelf  direc¬ 
tion  at  CHi2  were  properly  oriented  for  advecting  the  temp¬ 
erature  pattern.  These  cross-shelf  oscillations  were 
responsible  for  the  variability  in  the  bottom  temperature  4 
IP  after  the  passage  of  Frederic. 

C.  SOPEBI BESTIAL  RESPONSE 

is  Frederic  passed  within  130  km  to  the  west  of  the 
CHA1 ,  the  dominant  response  was  a  significant  increase  in 
the  westward  curmt  (Fig.  24).  This  increase  was  due  to 
convergence  of  flow  on  the  right  side  of  the  storm.  By  con¬ 
trast,  there  was  a  divergence  of  flow  on  the  left  side  of 
the  hurricane,  where  water  was  transported  off-shelf  by  the 
winds.  For  example,  a  negative  tide  was  recorded  at  Biloxi, 
Bississippi  which  was  indicative  of  the  left  side  of  the 


storm.  The  north-south  current  only  increased  by  about  20 
cb/s ,  whereas  the  east-vest  current  increased  to  80  cm/s. 
After  the  hurricane  passage,  near-inertial  waves  persisted 
over  the  subsequent  3  IP,  The  post-storm  semi-diurnal  tidal 
currents  were  acre  energetic  than  during  the  quiescent 
period  and  were  superposed  on  a  non-stationary  trend.  At 
the  intermediate  depth  at  CNA 1,  the  current  speed  increased 
to  80  cm/s  during  the  period  of  strong  forcing,  as  the 
semi-diurnal  tidal  currents  were  enhanced. 

The  presence  of  internal  tides  at  the  semi-diurnal  fre¬ 
quency  is  indicated  by  the  HKE  spectra  during  the  storm 
period  at  CHA 1  (Table  VIII) .  The  spectral  estimates  in  the 
semi-diurnal  frequency  band  are  almost  equal  to  the  inertial 
period  estimates.  Baines  (1973)  shows  that  one  necessary 
condition  for  the  presence  of  internal  tides  is  a  marked 
increase  in  the  semi-diurnal  tidal  currents.  The  observed 
increase  in  tie  semi-diurnal  tidal  currents  should  be  equal 
to  or  greater  than  the  inertial  effects.  The  storm-spectral 
estimates  at  the  semi-diurnal  frequency  exceed  those  from 
the  quiescent  period  by  acre  than  order  of  magnitude  (Appen¬ 
dix  B) .  A  barotropic  tide  propagating  over  the  shelf  break 
generates  internal  tides  cf  semi-diurnal  tidal  period 
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Figure  23. 


95  % 


cycles/hr 


Cross-Spectrue  Between  V-Velocity  and 
Teeperature  froe  c Ml 2  at  324  «.  K.  represents 
the  diurnal  tide  frequency,  £  the  inertial 
frequency  and  H.  the  seai-diurnal  tide 
frequency. 
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Figure  24.  Hid-depth  velocity  components  from  CHAl.  The 

ordinate  depicts  magnitude  in  cn/s  of  the  east- 
vest  currents  (upper)  and  north-south  currents 
(lover)  . 


(Prinsenberg  and  Rattray,  1975  ;  Barbee  §t  j^.,  1975  ;  Tor- 

\ 

grimson  and  Hickey,  1979).  Hence,  an  increased  barotropic  j 

tide  associated  with  the  storm  surge  propagating  over  rugged  ] 

bottom  topography  should  enhance  the  internal  tides  as  mani-  ; 

fested  in  the  semi-diurnal  tidal  currents.  * 
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TABLE  VIII 


Normalized  HKE  Spectra  at  Inertial  and  Semi-diurnal 

Frequency 


Heter 

Cepth 

(£) 


sail 


49 

64 

92* 


No.  of  Data 
Points 


Fr  equency 
He  solution 


(££h) 


Normalized 
HKE  Spectra 
f  M2 

(cm^s) 2/cp  h 


3200 

2410 

3000 


0.0018  2.0x10*3  1.0x10*3 
0.0025  3.0x10*3  6.0x10*2 
0  .0020  1.4x10*3  3.  2x10*3 


*  instrumentation  problems  with  time  clock 
H 2  is  the  semi-diurpal  tide  frequency  band 
f  is  the  inertial/diurnal  tide  frequency  band 


The  dceinance  of  internal  tides  on  the  shelf  at  chai  can 
be  explained  by  considering  the  bottom  slope  in  the  region. 
If  the  bottom  slope  is  greater  than  the  internal  wave  char¬ 
acteristic,  then  the  rays  will  not  be  admitted  to  the  shelf, 
but  instead  will  be  reflected  into  the  oceans  interior  (Bar¬ 
bee  gi  ai,  1975  ;  LeBlond  and  Hysaic,1978  ;  Torgrimson  and 
Hickey,  1979)  .  The  critical  slope  for  the  internal  wave 
characteristic  is  given  by: 


dz  ,  o 2  -  f 2 , 1  /2 
'  N2  -  a2 


(6) 


where  f  is  the  local  Coriclis  parameter,  which  is  equal  to 
7.18  x  10”*  rad/s  near  the  DeSotc  Canyon.  The  remainder  of 


the  variables  have  been  defined  above.  The  slope  of  the 
inertio-gravity  wave  characteristic  is  about  4  x  10-3, 
whereas  in  the  direction  of  the  canyon  axis  (45  deg.)  the 
bottoa  slope  is  6  x  10”3  .  Therefore,  it  is  clear  that  the 
inertio-gravity  waves  will  be  reflected  seaward  as  the  slope 
of  the  characteristic  is  less  than  the  bottom  slope  in  the 
DeSotc  Canyon  region.  Inertial  oscillations  will  be  locally 
generated  on  tke  shelf  but  they  will  only  persist  for  a  few 
IP  following  the  storm.  The  internal  tides,  however,  will 
be  admitted  to  the  shelf  region  because  the  slope  of  the 
semi-diurnal  tide  characteristic  is  2  x  10~2,  which  is 
almost  an  crder  of  magnitude  larger  than  the  bottoa  slope. 

The  semi-diurnal  tides  were  also  enhanced  at  CMA2  and  3, 
but  they  were  much  less  energetic  than  the  inertial  waves 
(Table  IX)  .  The  CW  polarized,  semi-diurnal  tidal  variabil¬ 
ity  was  more  than  an  crder  of  magnitude  greater  than  the  ccw 
variability  estimates,  except  in  the  bottoa  layers  at  CMA3. 
The  directions  of  the  semi-major  axes  of  these  waves  changed 
from  120  tc  20  deg.  with  depth  at  CHA2,  whereas  the  opposite 
occurred  at  CEA3,  where  the  direction  changed  from  30  to  50 
deg.  Hence,  the  waves  were  anisotropic  and  exhibited  a 
strong  north-south  component  of  flow  at  CHA2  and  to  the 
northeast  at  CBA3. 
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TABLE  IX 

Normalized  Hctary  spectra  at  the  Semi-diurnal  Frequency. 


Normalized 

Rotary 

Spectra 

Elli 

pse 

Depth 

CM 

c  cw 

Dir. 

Rotary 

(1) 

S)  2/CPh 

Stab. 

(Dea.) 

Coeff . 

a) 

ciy 

19 

1.3x10+3 

4.3x1 0+2 

0.80 

120 

♦  0.30 

179 

5.0x10+2 

9. 0x1 0+1 

0.73 

60 

♦  0.72 

324 

4.3x10+2 

1.2x1  0+2 

0.87 

20 

♦  0.41 

b) 

CBA3 

21 

1.0x10++ 

5.5x1  0+2 

0.58 

30 

♦  0.55 

251 

6.7x10+2 

1.0x10+2 

0.92 

90 

+  0.70 

437 

3.2x10+2 

1.  3x1  0+2 

0.40 

40 

♦  0.40 

457 

2.0x10+2 

1.  8x1  0+2 

0.15 

50 

♦  0.04 

In  summary. 

the  superinert ial 

response  in  the 

wake  of 

Frederic  was  dominated  by  the  enhanced  semi-diurnal  tidal 
current  at  all  three  arrays.  The  semi-diurnal  tidal  cur¬ 
rents  were  markedly  more  energetic  than  the  inertial  cur¬ 
rents  on  the  continental  shelf  at  CMA1 ,  and  agrees  with 
internal  tide  theory  (Baines,  1973)  .  The  storm  surge,  which 
is  associated  with  the  passage  cf  the  storm,  increases  the 
barotropic  tide.  As  this  increased  barotropic  tide  propa¬ 
gates  over  rugged  bottom  terrain  in  a  stratified  ocean,  the 
increase  in  tie  semi-diurnal  tidal  currents  are  manifesta¬ 
tions  of  internal  tides  (Earbee  et  jl,  1975  ;  Prinsenberg 
and  Battray,  1975  ;  Torgrimson  and  Bickey  ,  1979).  Internal 
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V.  NORMAL  BODES 


A.  THEORY 

The  Sturm-Licuville  problem  is  solved  using  the  vertical 
profile  of  Brunt-Vaisala  frequencies  and  linear  wave  theory. 
The  vertical  eigenvalues  cf  the  horizontal  velocity  are  then 
used  to  obtain  a  least  squares  fit  to  the  Fourier  coeffi¬ 
cient  time  series  computed  from  the  demodulation  of  the  cur¬ 
rent  meter  records.  The  amount  of  variability  associated 
vith  each  mode  is  determined  to  assist  in  understanding  the 
modal  structure  of  the  inertio- gravity  wave  response  gener¬ 
ated  by  the  passage  of  a  hurricane. 

The  problem  is  formulated  by  following  the  work  of 
Fjeldstad  (1958)  and  making  the  following  assumptions: 

•  f  plane; 

•  continuous  stratification; 

•  hydrostatic  balance  with  basic  state  at  rest; 

•  incompressible; 

•  inviscid;  and 

•  flat  bottcm. 
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The  governing  vertical  structure  equation  for  linear,  iner- 
tio-gravity  waves  is  given  by 

2 

^-4  +  U2$(Z)W  =  0  -)  (7) 

dz2 

where  H(z)  gives  the  vertical  structure  of  the  inertio-grav- 
ity  wave,  ,  *(z)*H2_  <r2,  and  K2  *  k2  ♦  l2,  and  k 

and  1  are  the  horizontal  wave  number  vectors  in  the  x  and  y 
directions,  respectively.  Furthernore,  equation  (7)  is  a 
second  order,  homogeneous,  ordinary  differential  equation. 
The  specification  of  the  problea  is  completed  by  imposing 
the  following  boundary  conditions: 

||  -  guW  *  0  z=0  7  (8) 

W  =  0  z=-d  •  (9) 

The  surface  boundary  condition  (8)  is  a  dynamic  boundary 
condition  which  states  that  at  the  sea  surface  pressure  is 
continuous  across  the  interface.  Boundary  condition  (9) 
specifies  nc  normal  flow  through  the  bottom.  The  solution 
to  equations  (7-9),  which  define  a  well  known  Sturm-Liou- 
ville  problea,  yields  the  structure  of  the  vertical  velocity 
for  constant  H,  and,  is  given  by 

Wn(z)  -  sin(^)  7  (10) 


where  the  mode  number  ns0,1,2....  The  horizontal  velocity 
eigenfunctions  0  (z)  are  given  b; 

Un(z)  ~  cos(^)  9  (11) 

which  is  the  vertical  derivative  of  (z)  . 

The  Sturm-Licuville  problem  is  solved  numerically  by  the 
predictor-corrector  fourth-order  method  modified  by  Hamming 
(Gerald,  1980).  This  method  computes  a  new  vector  from  four 
preceding  values,  A  fourth-order  Runge-Kutta  method  is  used 
for  the  adjustment  of  the  initial  vertical  increments  and 
the  computations  of  starting  values. 

The  resulting  amplitudes  of  the  horizontal  velocity 
eigenfunctions  fcr  modes  0  through  3  are  shown  in  Fig.  25  . 
The  number  of  zero  crossings  equals  the  mode  number.  Nest 
of  the  vertical  structure  lies  in  the  upper  thermocline 
region  near  the  base  of  the  mixed  layer. 

B.  DATA  ANAL1SIS 

The  horizontal  eigenfunctions  are  then  fitted  in  a  least 
squares  sense  to  the  time-varying  amplitudes  of  the  horizon¬ 
tal  velocities.  The  mathematical  details  are  given  in 
Appendix  A  with  the  final  matrix  equations  (16-18)  .  The 


time  evolution  of  the  north-south  velocity  coefficients  at 


Figure  25. 


Horizontal  Velocity  Eigenfunctions  at  29.33  N 
and  87.11  if.  The  abscissa  depicts  the  aeplitud 
of  the  horizontal  velocity  eigenfunction 
■ode  0,  (solid) .  lode  1  (dashed),  ** 

and  aode  3  (dashed-act ted)  . 


■ode  2 


_or 
(dotted 


CM 43  indicates  that  the  storm  excited  modes  0,1, and  3  with  a 
damped  mode  2  (Fig.  26)  .  The  amplitudes  of  modes  0  and  1 
reached  a  peak  on  JD  256  and  then  decreased.  By  JD  262,  the 
barotropic  mode  horizontal  speed  decreased  to  nearly  a  con¬ 
stant  value  of  about  10  cm/s.  After  JD  262,  modes  0  and  1 
oscillated  over  the  remainder  of  the  record  with  a  period 
of  about  5  IP.  This  behavior  of  the  modes  in  Fig.  26  indi¬ 
cates  that  a  modulation  occurred  between  inertial  wave  modes 
0  and  1.  More  importantly,  these  modes  define  a  modulation 
envelope  which  is  consistent  with  the  amplitude  variations 
of  the  currents  below  the  mixed  layer. 

At  CHA2,  the  time-varying  coefficients  of  the  modes  are 
markedly  different  in  their  behavior  (Fij.  27).  The  spatial 
variability  in  the  modal  structure  between  CM&2  and  C3A3  is 
associated  with  the  differences  in  the  mean  flow,  which  is 
topographically  controlled.  Initially,  all  of  the  modes 
were  excited,  with  a  fairly  energetic  mode  3.  The  baro¬ 
tropic  coefficient  decreased  to  about  2  cm/s  while  the  mode 
1  amplitude  approached  zero  near  JD  232.  After  JD  262, 
modes  0  and.  1  were  in  phase  and  defined  the  limits  of  the 


coefficients 


Figure  27.  Saee  as  Fig.  26,  except  at  CHA2 


The  medal  current  tine  series  are  re-coaputed  using  the 
coefficients  from  the  least  squares  fit  and  the  eigenfunc¬ 
tions  for  each  cf  the  current  aeter  depths  at  CM A3.  For 
example,  the  reconstructed  tine  series  are  compared  in  Fig. 
28  with  the  actual  demodulated  current  time  series  at  251  a. 
Recall  that  a  secondary  maximum  occurred  at  this  site  about 
10  IE  following  the  passage  of  Frederic  (Fig.  10) .  The 
reconstructed  amplitude  cf  the  barotropic  mode  is  greater 
than  the  actual  v  component  of  the  data  between  JD  255  and 
262  (Fig.  28).  However#  summing  modes  0,1  and  2  accounts 
for  most  cf  the  observed  variability#  with  the  barotropic 
mode  the  most  energetic.  Over  the  remainder  of  the  record, 
JD  262  tc  270#  most  of  the  observed  variability  can  be 
described  cnly  by  modes  0  and  2 .  The  contribution  of  the 
barotropic  mode  was  significant  and  consistent  with  earlier 
estimates  that  the  vertical  scales  were  of  the  order  of  the 
water  depth.  The  strength  of  the  barotropic  mode  continued 
for  the  entire  time  series  of  the  nodal  coefficients  at 
CMA3. 

Estimates  of  the  modal  contributions  to  the  observed 
near-inertial  variability  are  given  in  Table  I  for  both  the 
storm  and  poststorn  periods  at  CHA2#3.  These  estimates  are 
based  on  the  expression 


where  Vj  and  VjM  are  the  observed  and  nodal  scalar  conponents 
of  the  horizontal  velocity.  The  index  m  depicts  the  sunna- 
tion  of  the  nodes,  for  exanple,  0,  0+1,  and  0+1+2,  and  index 
j  repesents  a  sumnation  over  the  nunber  of  observations  in 
the  period.  The  storn  period  is  redefined  to  start  at  the 
tine  of  hurricane  passage  and  continue  for  7  IP.  The  post- 
storn  period  starts  where  the  storn  period  ends  and  contin¬ 
ues  for  roughly  6  IP  until  the  end  of  record. 

The  nixed  layer  variability  at  CH&2  is  doninated  by  the 
barotropic  node  during  the  storn  period.  The  contribution 
of  the  barotropic  node  tc  the  observed  horizontal  current 
variability  decreases  with  depth,  but  by  including  the 
baroclinic  nodes  most  of  the  near-inertial  variance  can  be 
described  by  nodes  0,1  and  2.  The  only  exception  is  at  179 
i.  The  first  three  nodes  contribute  only  54  and  37%  to  the 
observed  horizontal  current  variance  in  the  east-west  and 
north-south  directions,  respectively.  Furthernore,  the 


u-velocity  ccnpcnent  of  the  barotropic  node  exceeds  the 


post-storm  period  is  attributed  to  the  barotropic  mode.  The 
modes  are  such  more  well-behaved  during  this  period,  which 
is  presumably  a  manifestation  of  free  waves. 

The  modes  of  the  horizontal  velocity  are  well-behaved 
during  the  stcra  period  at  CMA3,  except  near  the  bottom 
where  only  20  to  22%  of  the  observed  variance  can  be  attrib¬ 
uted  to  the  barctropic  mode.  This  is  probably  due  to  bottom 
boundary  effects.  At  the  other  depths,  the  barotropic  mode 
explains  55  tc  73%  of  the  variance,  with  most  of  the  vari¬ 
ance  accounted  for  by  the  summation  of  the  barotropic  and 
first  two  barcclinic  modes.  During  the  post-storm  period, 
the  observations  of  the  mixed  layer  velocity  are  less  than 
that  cf  the  barctropic  mcde  in  the  north-south  direction, 
which  causes  the  value  to  be  less  than  zero.  The  reason  for 
this  large  discrepancy  is  due  to  the  e-folding  scale  of  the 
horizontal  velocity  in  the  mixed  layer,  the  persistence  of 
the  coefficient  cf  the  barotropic  mode  (see  Fig.  26) ,  and 
the  radiation  of  inertial  waves  away  from  the  storm  track. 
Otherwise,  73  tc  95%  cf  the  observed  horizontal  current  var¬ 
iance  can  be  attributed  to  the  barotropic  mode,  and  ;*dding 
the  first  two  baroclinic  mcdes  accounts  for  more  than  90%  of 


the  variance 


Hopkins  (1982)  modeled  only  the  baroclinic  response  zo 
hurricane  passage  and  he  found  that  the  inertio-gravity 
waves  decayed  too  rapidly  in  the  thermocline  and  bottom  lay¬ 
ers  to  account  for  the  vertical  structure  of  the  observed 
variability.  Thus,  a  model  with  only  the  baroclinic  modes 
could  not  explain  the  secondary  maximum  observed  in  these 
layers. 

Some  caution  has  to  applied  in  the  interpretation  of  the 
normal  modes  for  several  reasons.  First,  the  formulation  of 
the  Sturm-Liouville  problem  assumes  that  the  vertical  struc¬ 
ture  consists  of  a  standing  wave  for  the  N2  =  constant  case. 
From  basic  physics,  a  standing  wave  can  be  decomposed  into 
two  waves  of  equal  amplitudes  propagating  in  opposite  direc¬ 
tions,  which  assumes  no  phase  propagation,  over  the  first  7 
to  10  IP,  there  is  a  downward  propagation  of  energy,  as  well 
as  upward  propagation  of  phase,  which  violates  the  standing 
wave  argument.  Secondly,  a  flat  bottom  ocean  is  assumed  for 
the  application  of  normal  mode  theory.  The  validity  of  this 
assumption  crucially  depends  on  the  ratio  between  the  slope 
of  the  internal  wave  characteristic  to  the  bottom  slope. 
For  example,  as  this  ratio  approaches  unity,  the  energy  den¬ 
sities  in  the  bottom  boundary  layer  approaches  infinity.  Kt 
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Figure  28. 
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V. Component  Tiie  series  at  CH&3  for  the  observed 
time  series  (solid)  and  the  reconstructed  time 
series  (dotted)  for  a)  mode  0  ,  b)  mode  0+1  ,  c) 
mode  0+t+2  . 
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VI .  COHCLDSICNS 


The  ocean  response  near  the  DeSoto  Canyon  to  hurricane 
Frederic  vas  dominated  by  the  excitation  of  inertial  waves 
in  the  mixed  layer  and  inert io- gravity  waves  in  the  therno- 
cline.  These  waves  are  not  admitted  onto  the  shelf  because 


the  slope  cf  tie  characteristics  of  internal  wave  motion  is 
less  than  the  bottom  slope  (Barbee  et.  &1,  1975  :  LeBlond 

and  Bysak,  1S78  ;  Torgrimson  and  Bickey,  1979).  The  major 
features  of  the  forced  inertio -gravity  waves  observed  at 


CBI2  and  3  arc: 


i 


initially,  the 
the  order  of  1. 
the  water  column 
sage  cf  Frederic 


vertical  propagation  of  energy  was  of 
km/d  as  the  forcing  wa s  felt  throughout 
within  3  to  4.5  S  following  the  pas- 


•  the  initial  horizontal  propagation  of  energy  vas  abcut 
80  km/d  with  a  phase  speed  or  70  cm/s; 


•  the  waves  were  anisotropic  at  all  depths  and  had  an 
upward  propagaticn  of  phase; 


•  maximum  HRE  in  the  mixed  layer  occurred  at  the  fringe 
of  the  maximum  wind  regime  with  a  second  maximum  in  the 
thermocline  6  to  10  IE  after  the  storm; 


•  the  BKB  decayed  on  an  e-folding  scale  of  about  4  ip 
after  the  maximum  in  the  mixed  layer  and  the  second 
maximum  in  the  thermccline; 


•  the  vertical  propagation  after  about  6  IP  is  of  the 
order  of  10*  m/d  with  a  corresponding  horizontal  propa¬ 
gation  of  energy  of  30  km/d; 


•  the  nixed  layer  horizontal  scales  were  60  and  100  km  in 
the  cross  and  along-track  directions,  respectively, 

Ihile  the  scales  in  the  thermocline  were  25  and  50  km 
n  the  cross  ana  along-track  directions; 
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•  the  vertical  scales  were  cf  the  order  of  the  water 
depth  which  suggests  that  the  barotropic  node  was 
inportant; 


•  eost  of  the  observed  variabilty  was  in  the  barotropic 
■ode.  with  soae  contributions  froa  nodes  1  and  2; 

together  they  defined  a  aodulation  envelope  in  the  deep 
ther lccline  at  CBA3 ; 


•  the  nodal  coefficients  varied  in  space  as  well  as  tine 
due  tc  the  rugged  bottcn  topography  and  mean  currents. 


The  nean  flew  was  topcgrapically  controlled  and  reversed 
direction  fzen  surface  tc  botton.  The  forced  response  in 
the  nixed  layer  caused  the  nean  flow  to  increase  at  both 
CH k2  and  CHA3.  The  nean  flow  below  the  nixed  layer  was  also 
changed  by  the  store.  Furthermore,  this  nean  flow  was 
influenced  ty  sufcinertial  variability  (periods  of  3  IP)  and 
was  linked  to  the  advecticn  of  cooler  water  cross-shelf  rel¬ 
ative  to  CHI2. 

The  seai-diurnal  tidal  currents  increased  in  nagnitude 
and  varied  spatially  in  response  to  the  passage  of  hurricane 
Frederic.  The  near  dcainasce  of  the  seai-diurnal  tidal  cur¬ 
rents  over  the  inertial  actions  at  CHA1  was  consistent  with 
internal  tide  theories.  The  seai-diurnal  tides  were  adait- 
ted  tc  the  shelf  region  because  the  slope  of  the  internal 
wave  characteristic  exceeded  the  critical  bottoa  slcpe 


(Baines  ,  19*73  ;  Prinsenberg  and  Battray,  1975  ;  Torgrimson 
and  Hickey,  1979). 

In  a  broader  context,  the  oceanic  response  to  hurricane 
Frederic  was  a  geostrophic  adjustment  problem,  which  had 
both  transient  and  steady  state  components.  &fter  the  tran¬ 
sient,  inertic-gravity  waves  propagated  away  from  the  storm 
track,  a  gecstrophically  balanced  ridge  and  current  system 
remained  under  the  storm  track.  That  is,  the  steady  state 
currents  were  due  to  the  balance  between  the  horizontal 
pressure  gradients  and  the  Coriclis  force  (Seisler,  1970)  . 
Although  some  adjustments  to  the  mean  flow  are  described, 
additional  study  is  required  to  completely  understand  the 
interelaticnships  between  the  forced  mean  flow,  the  forced 
wave  field  and  the  bottom  topography. 


AFPEhDI  X  A 


k.  CCHPLII  01 BODULATIOH 

The  motivation  for  ccaplex  deaodulation  is  to  isolate 
the  carrier  wave  of  a  certain  frequency,  which  aust  be  known 
1  priori,  and  to  fora  amplitudes  and  phases  of  the  nodulated 
signal.  1  linear  filtering  approach  is  used  to  fora  the 
instantaneous  aaplitudes  and  phases  in  contrast  to  the  har- 
aonic  analysis  aethod.  Essentially,  the  anplitude  and  phase 
tiae  series  represents  a  local  haraonic  analysis  rather  than 
discrete  estiaates  averaged  over  continuous  periods  (Otnes 
and  Enochson,  1978)  . 

1  band-limited  tiae  series  is  aultiplied  by  the  trigono- 
aetric  arguacnts  of  the  carrier  wave.  The  resultant  tiae 
series  is  then  lew-pass  filtered  to  eliainate  soae  high  fre¬ 
quency  noise  that  is  generated.  The  cosine  and  sine  coeffi¬ 
cients  are  tfaer  coabined  tc  fora  the  aaplitude  and  phase  of 
the  nodulated  wave.  Consider  the  tiae  series: 

x  (i)*A  li)  cos  ( (2  v  ^i£t)  ♦  *(i))  ,  (1) 

where  4  is  the  carrier  frequency,  At  is  the  saapling  inter¬ 
val,  #  is  the  phase  ,  k  is  the  aaplitude  of  the  wave,  and  i 


represents  a  tiae  index.  Both  the  phase  and  asplitade  vary 
as  functions  cf  tiae.  Hultiply  equation  (1)  by  the  sine  and 
cosine  arguaents  of  the  carrier  frequency  fc  , 

X%  |i)*h  |i)|ccs  (( 2  ir  fe  iAt)  ♦  0  (i))^  sin  (2  ir  fc  iAt )  ,  (2) 

xc  |i)»h  (i)|ccs  ((2  IT  fe  iAt)  ♦  0  (i)  )j  cos  (2  n  fc  iAt )  ,  (3) 

where 

x,  (i)  *  x(i)  sin  (2irfciAt)  ,  and 
xe  (i)  3  x(i)  cos  (2irfeiAt). 

» 

Equation  (2)  can  be  rewritten  as 

x,  <i)  *  A  <i)/2  |  sin  <<47rfciAt)  ♦  #(i) )}-  sin  »  (i)  .  (4) 

Siailarly,  equation  (3)  can  be  expanded  and  simplified  as: 

xe  (i)  »  k  <i)/2  | cos  ( (4  ir  fc  iAt)  ♦  0  (!))}♦  cos  *(i)  .  (5) 

Low  pass  filtering  the  above  xs  ,xt  series  yields: 

y,  (i)  *  -1  (i)  /2  jsin  #  (i)J  ,  (6) 

ye  <i)  *  k  (i)  /2  | cos  *(i)|  , 

where  y4  <i)  and  ye(i)  are  the  low-pass  filtered  xs  (i)  and  x e 
(i) .  These  data  are  filtered  by  convolving  the  input  series 
with  a  set  cf  filter  weights  in  the  tiae  doaain  (Bendat  and 
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Piersol,  IS 7 1) .  Squaring  and  suaaing  ys  (i)  and  yc  (i) 
yields 

a*  li)  -  (y*  (i)  ♦  y |  <i>>  »  A*  /4 

or. 

Mi)  -  2  a  (i)  (7) 

where  Mi)  is  the  amplitude  as  a  function  of  tiae.  siai- 
larly,  the  phase  as  a  function  of  tiae  is  coaputed  as: 

♦c(i)  »  Tan-»  <-ys  (i)  /  ye  (D )  .  (8) 

B.  E0TBNTI1L  ENEBGT 

The  potential  energy  is  coaputed  to  deteraine  the  rela¬ 
tionship  of  the  isothera  displaceaents  to  internal  wave 
aoticn.  Given  a  teaperature  tiae  series  froa  a  aoored  array 
and  the  corresponding  vertical  teaperature  gradient  df/dz, 
the  isothera  displaceaent  t  is  coaputed  in  the  following 
Banner: 

c.  .  T—  ,  O) 

dT/dz 

where  T'  is  the  fluctuating  part  of  the  observed  teapera- 
tures.  The  potential  energy  per  unit  voluae  is  defined  by: 


(10) 


PE  » 


P°  M2 
2~  N 


share  H*  is  the  firunt-Vaisala  frequency#  f*  is  the  isothera 
displaceaent  fluctuations#  and  p.  is  the  reference  density, 
which  is  assused  to  be  unity. 

C.  LEAST  SQU1BES  FIT 


Given  a  demodulated  tine  series  of  the  aaplitude  coeffi¬ 
cients  and  the  vertical  structure  of  the  horizontal  velocity 
eigenfunctions#  the  probles  is  to  cospute  the  time-varying 
coefficients  associated  with  each  dynamical  mode.  The 
time-varying  amplitudes  0(2^#t)  at  depths  z,,z2#  ....  z* 

may  be  written  as 


— 

0  (zi  ' t)m 


where  EM(Z|)  are  the  eigenfunctions  for  modes  1,2,3.  ..#  and 
C„<t)  are  the  coefficients  to  be  calculated  in  the  analysis. 
Only  three  fcarodinic  medal  coefficients  are  computed 
because  current  meters  were  deployed  at  three  depths  in  the 
water  column.  Por  example#  the  current  at  depth  z,  ,  is 


written  as 


0,  *  C,  E,  <Z,  )  ♦  C,Ea(Z,)  ♦  CjEj  (Z|  )  # 
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and  theca  are  siiilac  expressions  at  depth  z2  and  z3.  The 
square  of  the  error  in  the  least  squares  fit  at  depth  z,  is, 

e*  *  |o,  -  (C(  E|  (Z|  )  ♦  C» B* (Z| )  ♦  CsE3(z,))}*  (11a) 

and  sinilarly  at  z  and  z  , 

e|  *  |oa-  (C,  E,  (za)  ♦  CaEa(za)  ♦  CaEa(za))|*  ,  (11b) 

ef  *  {03-  (C,  E|  (za )  ♦  CaE4(z3)  ♦  C1E>(zJ))J*  .  (11c) 

The  error  is  einisized  by  solving  the  following  set 

def  /*CX  *de*  /dC,  +  de|  /dc,»  0  ,  (12a) 

def  /dCa  ♦de*  /dCa+de|  /dcx»  0  ,  (12b) 

def  /c)C3  ♦de/  /dC3+de|  /dCs*  0  .  (12c) 

Differentiating  equation  (11)  with  respect  to  c,  and  sunning 
the  expressions  according  to  equation  (T2a)  yields 


0,  I,(z,)  ♦  OaEj(za)  ♦  D3Ej(z3)  » 

C,(E,(Z,)  E  |  (Z ,  )  ♦El(Zl)  E,  (Z2)  *  E  (  (z3 )  E,(Zj)) 


♦  Cj  (E(  (z, )  E a  (Z  |  )  ♦  E  |  (za  )  E  j  (Za)  ♦  E,  (Zj )  Ej  (Za)  ) 

♦  Cj  (E,  (z, )  E3(Z,)  ♦  E,(za)  E3(za)  ♦  E  |  (zs)  EJ(z3)), 


or 


£o{I,  (*i)  »  C,£e,  (Zi)  E  |  (Z  ( )  ♦CaVE,(zl)  Ex(zi 

Tl  tM  Ul 
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♦  <Z i  >  ES<*i>  •  03) 

Mow  if  equations  (12b,c)  are  applied  to  equations  (11)  ,  the 
corresponding  equations  for  C2  and  Ca  become: 

*  C.^B,  (2,1  !*<*•,)  ♦  ca^Ea(zi  >  Ea  (zi  > 

♦C^S^z,)  Ea(Zf)  ,  (14) 

^SjjIa(Zj)  -C^E.fZ;)  E,<Z.)  *  C2^E2(Zi)  By  (Z  j  ) 


*  c3yj2(z\)  E a ( z i ) 


By  inspection,  let 


<Wa^EM<2i 


)  Em(Z() 


for  a,n  =1,2,3  ,  (16) 


*.*<*;> 


Therefore,  equations  (13-15)  can  be  simplified  using  the 
above  definitions  into  a  simple  matrix  equation  of  the  fora 


&+  *  aMM  * 


MU 


Equation  (18)  is  a  simple  symmetric  matrix  equation  and  it 
was  solved  using  the  IHSL  routine  LEQT2F  on  the  IBM  3033. 


APPENDIX  E 


A.  DEEP  BZGICB 

The  ocean  currants  at  CHA2,3  during  the  quiescent  period 
clearly  showed  topographical  effects  associated  with  the 
OeSotc  Canyon.  The  behavicr  of  the  surface  currents  is  best 
described  by  a  FVD  (Fig.  29).  The  strong  u  coaponent  at 
CHA2  produced  a  total  displacement  of  600  ka  over  39  days 
for  a  15.4  ka/d  mean  current.  Inertial  waves  were  super¬ 
posed  on  this  aean  flow,  particularly  during  the  period  froa 
JO  232  to  239.  The  aean  flow  increased  during  this  period. 
At  CHA3,  the  aean  flow  decreased  in  the  southerly  direction 
froa  13.3  to  9.5  ka/d  (Fig.  29b)  .  The  different  orientation 
in  the  aean  flow  across  the  entire  shelf  region  is  presum¬ 
ably  due  to  the  influence  of  bottoa  topography. 

The  trajectory  of  the  near-bottoa  horizontal  currents  at 
CHA2  (Fig.  30a)  was  in  the  saae  direction  as  in  the  nixed 
layer  until  JD  232.  After  the  reversal,  the  aean  flow  was 
towards  the  west  at  about  3.3  ka/d.  The  bottoa  PVD  at  CHA3 
(Fig.  301)  showed  that  initially  the  aean  flow  was  towards 
the  north  at  2.0  ka/d  over  the  first  6  IP.  The  flow  then 
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Figure  29. 


23  ®T  a 

Surface  pvd  fcr  a)  C HA2  and  b)  CHA3.  The  * 
represents  a  tine  span  of  48  h  starting  a)  JD 
2C9  for  CHA2  (19  m)  and  b)  JD  213  for  CHA3  (21 
»)  - 


>2 


U’4.1 


changed  direction  and  accelerated  tc  about  2.8  km/d  in  the 
saae  direction  as  in  the  surface  layer. 

Horizontal  kinetic  energy  spectra  for  this  period  are 
summarized  in  Table  XI  for  both  inertial  and  semi-diurnal 
frequency  hands.  The  vertical  variation  in  the  inertial  HKE 


levels  at  CHA2  indicated  a  larger  amount  of  energy  in  the 
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Pigute  30.  Bcttcm  PVD  for  a)  CH A2  and  b)  CMA3.  The  * 

represents  a  tine  span  of  48  h  starting  a)  JC 
209  for  CHA2  <324  a f  and  b)  JD  213  for  CHA3  (457 
■)  . 

theraocline  than  at  the  surface  or  near-bottom.  However, 
the  difference  between  the  surface  and  bottom  energy  levels 
was  about  cne  order  of  magnitude.  At  CHA3 ,  the  inertial  HKE 
estimates  decreased  steadily  with  depth  and  there  was  more 
than  an  order  of  magnitude  difference  between  the  surface 
and  bottom  layers.  The  HKE  levels  were  slighty  greater  at 


.A] 


ir  * 

<: 

m 


CH  A3  in  tee  surface  layer  and  theraocline,  but  they  were 
slighty  less  near  the  bottca  than  at  CBA2. 


TABLE  XI 

Horaalized  HKE  Spectral  Estiaates  at  9  5%  Confidence 


Beter 

Depth 

<1) 


No.  of  Data 
Points 


Bandwidth 


Noraalized 
HKE  Spectra 
f  B2 

(ca/s)  ♦2/cph 


a 

?  W 

cm 

% 

« 

« 

19 

5780 

0.001 

1.7x10+3 

2. 1x10*2 

4 

* 

179 

"590 

0.0013 

3. 0x10*3 

3.8x10*1* 

3 

1 

J 

324 

5780 

0.001 

3. IxlO** 

8. 0x10*1* 

§ 

!  fe) 

) 

CJS13 

j 

»  » 

21 

4680 

0.0013 

8. 0x10*3 

2.3x10+3 

j 

251 

4680 

0.0013 

4.5x10*3 

7.0x10+2 

*53 

1 

i 

437 

4680 

0.0013 

1.9x10*3 

4. 6x10+3 

457 

4680 

0.0013 

2.0x10*2 

1.7x10+1 

*  instruaentaticn  prebleas  with  tiae  clock 
H2  is  the  seei-diurnal  tide  freguency  band 
f  is  the  ineitial/diurnal  tide  frequency  band 


The  seai-diurnal  HKE  estiaates  also  varied  spatially. 
The  estiaates  cf  HKE  in  the  theraocline  and  near-bottoa  lay¬ 
ers  at  CHA2  were  not  significant;  however#  the  HKE  estiaates 
at  CHA3  were  significant  at  all  depths.  The  HKE  estiaates 
at  437  a  even  exceeded  both  the  surface  layer  and  therao¬ 
cline  levels  by  about  an  order  of  aagnitude.  In  coaparison 
with  the  near-bottoa  estiaates  at  457  a#  the  kinetic  energy 


levels  were  acre  than  two  orders  of  nagnitude  greater. 
These  pronounced  changes  in  the  seni-diurnal  energy  levels 
with  depth  indicate  the  presence  of  internal  tide  notion. 
In  aost  of  tie  kinetic  energy  spectra  exanined,  longer 
pericd  notions  were  insignificant,  because  their  energies 
were  saeared  acrcss  adjacent  bands  except  in  the  bottoa  lay¬ 
ers.  Subinertial  variability  at  periods  of  about  3  IP  was 
quite  apparent  in  the  bottca  records  at  both  CMA2  and  CH&3. 
These  notions  were  generally  nore  energetic  than  the  seni- 
diurnal  tides. 

notary  spectra  are  analyzed  for  the  surface  and  botton 
current  records  and  are  given  in  Table  XII.  Host  of  the  hke 
variability  is  consistent  with  polarized,  Cl  rotating  iner¬ 
tial  notions.  The  Cl  energy  level  exceeds  the  CCS  level  by 
about  an  order  of  nagnitude.  The  ellipse  orientations  for 
these  waves  rotates  CCS  (hacks)  with  depth.  These  ellipses 
are  stable  except  for  the  near-botton  inertial  period  notion 
at  CBA2.  The  direction  of  the  seni-najor  axis  of  the 
ellipse  varies  over  tine,  which  indicates  that  the  inertial 
waves  are  isotropic.  The  seni-diurnal  tidal  notion  is 
polarized  in  the  Cl  direction.  Generally,  the  seni-diurnal 
tidal  notions  are  coherent,  although  the  phase  differences 
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bet v ten  the  horizontal  velocity  coaponents  vary  between  the 


arrays. 


TIBLE  XII 

Serialized  Rotary  spectrum  Analysis 


Free.  Depth 

llAi  (a) 


Iforaalized 
Rotary 
Spectra 
C9  CCW 

(SI^S)  tiissk 


Ellipse 

Stab.  Dxr.  Rotary 
CSSS*)  ££&££• 


a)  5JSI2 


I 

19 

1.6x10** 

1.9x10** 

0.80 

103 

♦0.90 

s 

19 

1.6x10** 

3.6x10*  * 

0.39 

101 

♦0.87 

I 

324 

5.0x10** 

1.4x10*  * 

0.05 

90 

♦0.87 

s 

•  324 

1.0x10*0 

1.1x10*o 

0.60 

169 

-0.18 

k)  £MA2 


i 

s 

i 

s 


1.  2x1 0** 
4. 0x1  0** 

hmtt 


1.7x10** 
3. Ox  10*  1 
1.6x10** 
3.9x10*® 


0.40 

0.61 

0.51 

0.20 


157 

55 

T 


♦0.99 

♦0.63 

JM! 


I»inertial  frequency  band 
S*Seei-diurnal  frequency  band 


Ocean  bottca  temperatures  ranged  between  8.0  and  12.5  °c 
at  CBA2  and  8.0  to  9.0  ®C  at  CBA 3  (Fig.  31) ,  but  there  were 
different  variations  and  trends.  At  CBA2,  (not  shown)  the 
teaperature  decreased  gradually  to  about  8.3  oc  near  JD  234, 
and  then  increased  abruptly  without  any  periodicities.  Bct- 
tos  temperatures  at  CBA3  were  such  lover  and  there  was 
near-inertial  aotion  superposed  on  this  early  cooling  trend. 
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The  ainiaua  teeperature  occurred  on  JO  232,  which  was  about 
2  IP  before  tbe  ainiaua  observed  at  CHA2.  These  ainiaa  in 
the  teaperature,  in  addition  to  the  cooling  trends,  are  man¬ 
ifestations  of  the  advection  of  cold  water  froa  the  DeSotc 
Canycn  by  tbe  aean  flew. 

B.  SBALLOB  BEGICN 

In  the  theraocline  (49  a) ,  tbe  inertial  oscillations 
were  not  nearly  as  energetic  and  did  not  rotate  cw  as  in  the 
aired  layer  (Fig.  32)  .  The  non-stationary  trends  in  the 
east-west  currents  were  siailar  to  the  aired  layer  observa¬ 
tions,  but  there  was  an  initial  aean  flow  towards  the  deep 
ocean.  Superposed  on  the  aean  flow  were  lower  frequency 
oscillations  with  periods  cf  the  order  of  about  10  to  14  IP. 
The  total  displaceaent  over  the  forty  days  was  about  300  ha 
which  corresponded  to  a  aean  current  of  7.5  ka/d. 

Inertial  oscillations  at  64  a  were  not  as  energetic  as 
at  49  a,  but  there  was  a  siailar  non-stationary  trend  in  the 
tine  series.  This  trend  was  associated  with  a  lower  fre¬ 
quency  waves  having  a  period  of  10  to  14  IP.  These  oscilla¬ 
tions  were  superposed  on  a  northward  aean  flow  which  was 
sligbly  less  than  7.0  ka/d.  Records  froa  the  current  aster 
nearest  to  the  bottoa  (9  2  a)  ,  were  out  of  synchronization 


Figure  31.  Teeperature  Tiee  Series  froe  CMA3  at  457 


Figure  32 


Bid-depth  Pft  at  CH11.  The  *  represents  a  tiee 
span  of  4  8  h  starting  on  JD  207  at  49  ■. 


doe  to  internal  clock  probleas.  For  exaaple,  large 
increases  in  the  velocity  occurred  on  JD  252  as  opposed  to 
256.  Hence#  the  data  were  not  included. 

The  HKE  spectra  for  the  21#  49#  and  64  a  depths  are 
shown  in  Table  XHX  for  both  the  inertial  and  seai-diurnal 
frequency  bands.  The  inertial/ diurnal  tidal  kinetic  energy 


in  the  surface  layer  was  two  orders  of  aagnitude  larger  than 
in  the  subsurface  layers.  The  semi-diurnal  tidal  motion 
decreased  from  the  surface  to  49  ■  and  then  increased  with 
depth  such  that  the  HKE  near  the  bottom  almost  equaled  that 
near  the  surface. 


tarts  decreased  to  a  relative  niniaun  as  cold  vater  was 
advected  by  tbs  scan  flow. 

In  suaaary,  based  on  the  observations  daring  the  quies¬ 
cent  period,  the  aean  flew  strcngly  depends  on  the  bottoa 
topography.  Inertial  band  notion  is  evident  at  all  depths 
daring  the  quiescent  period.  Scae  of  this  notion  is  due  to 
freely  prepagating  inertial  waves  while  the  renainder  is 
part  of  the  forced  diurnal  tides.  Internal  tides,  as  well 
as  inertial  oscillations,  are  part  of  the  freely  propagating 
internal  wave  continuua  below  the  surface  layer.  The  effect 
of  the  topography  on  the  aean  flow  and  the  wave  notions  fur¬ 
ther  coaplicates  the  dynaaics  of  the  region  and  the  resul¬ 
tant  circulation  .  Longer  period  waves  (2.5  to  3  day 
period)  set  up  north-south  oscillations  that  advect  cooler 
water  froa  the  Desoto  Canyon  onto  the  shelf  region.  Super- 
inertial  frequency  notions  are  influenced  by  the  diurnal  and 
seai-diurnal  tides  and  contribute  to  the  variability. 
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